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Advanced  Research  Workshop  on  Quantum 
Transport  in  Semiconductors 

Hotel  Villa  del  Mare,  Maratea,  Italy  17-22  June  2001 


Session  1 

Moderator:  John  Barker,  Glasgow  University,  "Status  of  the  physics  and  modeling  of  ultrasmall 

devices" 

•  1 A  Greg  Timp,  University  of  Illinois,  "Making  small  MOSFETs:  ballistic  and  quantum  effects" 

•  I  B  Max  Fischetti,  IBM  Research,  "Modeling  small  MOSFETs-the  role  of  quantum  and 
many-body  effects" 

•  1C  Steve  Goodnick,  Arizona  State  University,  "Full  Band  Stmcture  Calculations  for  Transport 
in  Wide-Band-Gap  Semiconductors" 

Session  2 

Moderator:  Max  Fischetti,  IBM  Research,  "Semi-classical  modeling  of  small  semiconductor  devices" 

•  2A  Asen  Asenov,  Glasgow  University,  "Discrete  impurities  and  quantum  potentials  in 
MOSFET  modeling" 

•  2B  Richard  Akis,  Arizona  State  University,  "Effective  potentials  for  quantum  effects  in 
MOSFETs" 

Session  3 

Moderator:  Antti-Pekka  Jauho,  Danish  University  of  Technology,  "Introduction  to  Quantum 

Transport" 

•  3A  David  Feny,  Arizona  State  University,  "Wave  function  approaches  for  self-consistent 
computations  of  transport  in  quantum  dots  and  arrays" 

•  3B  John  Barker,  Glasgow  University,  "Trajectories  in  quantum  mechanics" 

•  3C  Michael  Bonitz,  Rostock  University,  "Non-equilibrium  Green's  functions:  transient 
phenomena  and  the  role  of  the  initial  state  for  devices" 

Session  4 

Moderator:  Chihiro  Hamaguchi,  Osaka  University,  "The  Metal-Insulator  Transition" 

•  4A  Gunther  Bauer,  University  of  Linz,  "The  metal-insulator  transition  in  d  =  2" 

•  4B  Jonathan  Bird,  Arizona  State  University,  "The  metal-insulator  transistion  in  open  quantum 
dots  and  arrays" 

Session  5 
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Moderator:  Michael  Bonitz,  Rostock  University,  "Spins  and  few  particle  problems  in 
Semiconductors" 

•  5A  David  Awschalom,  UC  Santa  Barbara,  "Spin  coherence  and  optical  measurements" 

•  5B  Daniel  Loss,  Uni  versity  of  Basel,  "Quantum  computing  in  semiconductor  systems" 

•  5C  Sankar  Das  Saima,  University  of  Maryland,  "Few  elecUon  (and  few  impurity)  systems-small 
system  effects" 

Session  6 

Moderator:  Gerhard  Klimeck,  Jet  Propulsion  Laboratory,  "Applications  of  quantum  transport  in 
devices" 

•  6 A  Carlo  Jacoboni,  Modena  University,  "The  Wigner  function  and  quantum  transpoiT' 

•  6B  Harold  Grubin,  SRA,  Inc.,  "Modeling  resonant  tunneling  diodes  with  Wigner  functions  and 
density  matrices" 

•  6C  Dejan  Jovanovic,  Motorola,  "Non-equilibrium  Green's  functions  for  MOSFET  modeling" 

Session  7 

Moderator:  Karl  Hess,  University  of  Illinois,  "Applications  of  quantum  transport  in  optics" 

•  7A  Tilmann  Kuhn,  University  of  Munster,  "Quantum  kinetics  and  the  femtosecond  time  scale  in 
optical  excitation  of  semiconductors" 

•  7B  Rolf  Haug,  University  of  Hannover,  "Single-electron  charging  effects  in  quantum  dot 
arrays" 

Session  8 


Moderator:  Anant  Anantram,  NASA  Ames,  "Novel  New  Concepts" 


•  8  A  Karl  Hess,  University  of  Illinois,  "EPR  Experiments:  the  Assumptions,  and  the  Failure  of 
Bell's  Theorem" 

•  8B  Antti-Pekka  Jauho,  Danish  University  of  Technology,  "Scanning  probe  measurements  of 
biomolecules  on  silicon  surfaces" 

•  8C  Mark  Lundstrom,  Purdue,  "The  Landauer  approach  in  device  modeling" 

•  8D  Wolfgang  Windl,  Motorola,  "Diffusion  and  clustering  of  impurities-a  problem  that  cannot 
be  ignored" 

Poster  Presentations 

Poster  Group  1 

•  1-1  "Analytical  theory  for  the  low  frequency  transport  and  noise  in  QID  conductors,"  G. 
Gomilla  and  L.  Reggiani,  Universidad  Polytechnica  Catalunya. 

•  1-2  "Monte  Carlo  simulations  of  Wigner  function  tumieling:  Role  of  the  effective  potential,"  L. 
Shifren  and  D.  K.  Feny,  Arizona  State  University. 

•  1-3  "3D  device  Monte  Carlo  modeling  with  discrete  impurities,"  S.  Bannud,  S.  Galdin,  P. 
Dollfus,  Institut  d'Electronique  Fondamentale. 

•  1-4  "3D  modeling  of  imper  fect  interfaces  and  edges  in  decanano  MOSFETs,"  S.  Kaya,  Ohio 
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University,  and  A.  R.  Brown,  S.  Roy,  and  A.  Asenov,  Glasgow  University. 

•  1-5  "Comparison  of  first-order  quantimi  correction  schemes  in  3D  drift-diffiision  simulations  in 
sub-0.1  micron  MOSFETs,"  J.  R.  Watling,  A.  R.  Brown,  A.  Asenov,  and  D.  K.  Ferry,  Glasgow 
University. 

•  1-6  "Quantum  potential  in  the  presence  of  heterointerfaces,"  .1.  R.  Watling,  R.  C.  W.  Wilkins, 
and  J.  R.  Barker,  Glasgow  University. 

•  1-7  "Performance  of  scaled  double  delta  doped  PHEMTs,"  K.  Kalna,  J.  R.  Watling,  A.  Brown, 
and  A.  Asenov,  Glasgow  University. 

•  1-8  "Electron  transport  and  infrared  photoconductivity  in  quantum  dot  structures,"  V.  Ryzhii 
and  V.  Mitin,  University  of  Aizu. 

•  1-9  "Photoluminescence  from  hot  electrons  in  low  dimensional  systems,"  H.  Moraose,  Y.  Inui, 
M.  Itoh,  and  C.  Hamaguchi,  Osaka  University. 

•  1-10  "Nonlinear  transport  through  an  array  of  quantum  dots,"  G.  Kiesslich,  A.  Wacker,  and  E. 
Scholl,  Technische  Universitat  Berlin. 

Poster  Group  2 

•  2-1  "Quantum  algorithms  in  the  frame  of  a  coupled-quantum- wires  physical  system,"  S. 
Reggiani,  A.  Bertoni,  R.  Brunetti,  and  M.  Rudan,  University  of  Bologna. 

•  2-2  "Quantum  waves  in  the  coupled-wire  Qubit,"  J.  Elarris,  R.  Akis,  and  D.  K.  Ferry,  Arizona 
State  University. 

•  2-3  "Quantum  mechanical  study  of  nanoscale  MOSFET,"  A.  Svizhenko,  M.  P.  Anantram,  and 
T.  R.  Govindam,  NASA  Ames  Research  Laboratory. 

•  2-4  "Effects  of  lead  population  on  magnetic  flux  controlled  dissipative  electron  transport 
through  coupled  quantum  dots,"  N.  Horing,  Stevens  Institute  of  Technology. 

•  2-5  "Transport  in  quantum  dot  arrays,"  N.  Mori,  Y.  Takamura,  T.  Ishida,  and  C.  Hamaguchi, 
Osaka  University. 

•  2-6  "Non-equilibrium  transport  in  nanostructures,"  J.  Fransson,  University  of  Uppsala. 

•  2-7  "Tunneling  spectroscopy  of  exchange-correlation  interaction,  of  electrons  in  Schottky  barrier 
in  quantizing  magnetic  field,"  A.  Shulman,  Institute  of  Radioengineering  and  Electronics. 

•  2-8  "Stochastic  simulation  of  the  Barker-Ferry  equation,"  M.  Nedjalkov,  Teclinical  University 
of  Vienna. 

•  2-9  "Simulation  of  entanglement  in  semiconductor  quantum  wires,"  A.  Bertoni,  R.  lonicioiu,  P. 
Zanardi,  F.  Rossi,  and  C.  Jacoboni,  University  of  Modena. 

•  2-10  "Mesoscopic  fluctuations  of  Coulomb  drag  between  quasi-ballistic  I  D  wires,"  N.  A. 
Mortensen,  K.  Flensburg,  and  A.-P.  Jauho,  Danish  Technical  University. 

Poster  Group  3 

•  3-1  "Spin  dynamics  in  IIl-V  quantum  wells  using  Monte  Carlo  simulation,"  A.  Bournel  and  P. 
Hesto,  Universite  Paris  Sud. 

•  3-2  "Electron-spin-phonon  coupling  and  relaxation  dynamics  in  a  double  quantum  dot,"  V. 
Puller,  L.  Mourokli,  and  N.  Horing,  Stevens  Institute  of  Technology. 

•  3-3  "Quantum  computing  with  spin  qubits  in  semiconductors  structures,"  V.  Privman,  Clarkson 
University. 

•  3-4  "High-field  transport  thi'ough  semiconductor  heterostructures,"  M.  Morifuji,  Osaka 
University. 

•  3-5  "Effect  of  electronic  disorder  on  the  phonon-drag,"  V.  Mitin,  Wayne  State  University. 

•  3-6  "Wi.gner  function  dynamics  in  presence  of  an  infinite  potential  barrier,"  P.  Bordone  and  C. 
Jacoboni,  Modena  University. 

•  3-7  "Quantum  transport  in  2D  MOSFETs,"  A.  P.  Anantram,  A.  Svizhenko,  and  T.  R.  Govindan, 
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NASA  Arnes  Research  Center. 

•  3-8  "Density-matrix  modeling  of  terahertz  photon-assisted  tunneling  in  resonant  tunneling 
diodes,"  M.  Asada,  Tokyo  Insitute  of  Teclmology. 

•  3-9  "Hyperfme  interaction  between  nuclear  and  electronic  spins,"  I.  D.  Vagner,  Holon 
Academic  Institute  of  Technology. 

•  3-10  "Metal -insulator  transition  in  2D  few-electron  systems,"  A.  Filinov,  M.  Bonitz,  and  Yu. 
Lozovik,  University  of  Rostock. 

•  3-11  "Quantum  dot  modeling  using  NEMO-3D,"  G.  Klimeck,  F.  Oyafuso,  R.  C.  Bowen,  T.  B. 
Boykin,  T.  A.  Cwik,  E.  Huang,  and  E.  Vinyard,  Jet  Propulsion  Laboratory. 
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ultra-small 

Semiconductor  devices 
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The  dinosaur  model:  silicon  will  eventually  become  extinct 


Time  (challenged  by  JRB/DKF 
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(International  Roadmap  for  Semiconductors  - 1999  Edition) 


•  What  are  the  implications  of  the  dinosaur  refusin 
to  die  for  quantum  devices? 


Certainly,  below  50  nm  we  expect  strong 

atomistic  effects. 

The  predicted  high  carrier  densities  also 
suggest  strong  many-body  effects:  Coulombic, 
correlations/fluctuations. 
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The  Secret  Behind  Moore’s  Law 
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2000  ITRS  National  Roadmap  Projections 


30nm  research  is  targeted  at  a  ~15  year  horizon 

Technology  Generations 


YEAR  OF  FIRST  PRODUCT  SHIPMENT 

2001 

^^9  2008 

2014; 

TECHNOLOGY  GENERATIONS 

DENSE  LINES  (DRAM  Half-Pitch)  (nm) 

130 

90  Igl 

30 

Logic  (HIgh-Volume:  MPU) 

Logic  Transistors/chip  (including  SRAM) 

67.3M 

■Baal  539M 

4;3B 

Chip  Frequency  (GHz) 

On>chip  local  clock,  high  performance 
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TECHNOLOGY  REQUIREMENTS 
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160 
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*Drive  current  independent  ot  technology 
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•exponential  increase  in  leakage  with  increasing  Vq 
•exponential  increase  in  leakage  with  tg^ 


OXIDE  THICKNESS(nm) 


Specimen  thickness 


•  Apparently  very  sharp  interface  !I 

Rough  interface  (11  A  pp)  looks  as 
smooth  as  atomically  sharp  interface. 
Minor  Intensity  differences 
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Scanning  Transmission  Electron  Microscopy/ 
Electron  Energy  Loss  Spectroscopy 

X-Ray  Detector 
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EELS  and 
Bright  Field 
Detectors 

Increasing 
Energy  Loss 


What  does  t^^=1nm  look  like? 
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Nominally  1.1  nm  thick  gate  oxide:  0.8-1.0nm  Si02 


1 .6nm  wide  oxygen  profile 
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The  MOSFET  from  Landauer’s  Perspective 
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EFFECTIVE  MOBILITY,  (cmWs) 


MOBILITY  RESOLVED 
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S.  Goodnick 


Roughness  limited  mobility  decreases  as  inverse  square 
of  surface  field,  and  as  (AA)^: 

>  Si  100  Roughness  Limited  Mobility 


S.  Goodnick 


Scanning  Tunneling  Microscopy 
THROUGH  the  Oxide 


ELECTRON  MOBILITY  (cmWsec) 


AA=0.67nm-2 1=>  0.35nnT2 


K.  Evans-Lutterodt 


MOBILITY  MEASUREMENTS  VS.  CALCULATION 


A=0.28nm,  A=2.4nm 


Completely  constrained  mobility  Calculation  on“rough 
MOSFETS  correlates  well  with  measurements. 


•  Homeric  challenges  face  bulk  silicon  technology,  but  the 
same  is  true  of  all  the  pretenders  to  the  throne. 

•  Some  of  these  challenges  stem  from  “atomic”  sensitivity 
(e.g.  the  gate  oxide) 

•  As  process  control  decreases;  circuit  margins  must 
increase  at  an  even  faster  rate. 


Long-range  Coulomb  interactions, 
coupled  plasmon /insulator-T 0-phonon 
modes,  and  electron  mobility  in  Si 
inversion  layers 


Max  Fischetti 


IBM  Research  Division 
T.  J.  Watson  Research  Center 
Yorktown  Heights,  NY  10598 

in  collaboration  with  Deborah  Neumayer  and  Eduard  Cartier 


June  2001 
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Soft  bonds  imply  ‘soft’  TO-modes.  huJro  ~  energy  low  enough 

to  allow  emission  by  thermal  electrons,  low  enough  to  be  thermally 
excited  (absorption  possible). 
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Effective”  dielectric  functions  RPA  vs  long-wavelength-limit  for  Si/Si02/Si 
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Coupled  modes  for  Si/Si02/SI  system:  Coupled  modes  for  Si/Si02/Si  system: 

Effect  of  coupling  Thickness  dependence 
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Coupled  modes  for  Si/AbOs/Si  system  Coupled  LO-plasmons  modes  for  Si/AIN/Si 

system 
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Coupled  LO-plasmons  modes  for  Si/ZrOa/Si  TO-  and  plasmon  content 

system  Landau  damping 
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FTIR  spectra  (D.  Neumayer)  Mobility 


WAVENUMBER  (cm"'')  "matrix  form: 


Scattering  processes  •  Surface  roughness: 

Ando's  model  assumes  Ist-order  pertubation  given  by  steps  at  interface  of 
rms  height  A  and  correlation-distance  A. 
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Subbands  Triangular-well  approximation 
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The  Si/Zr02/Si  system  Effective  mobility  I  (average  Ug) 


■PLANE  WAVE  VECTOR  Q  (10®  cm“^)  IN-PLANE  WAVEVECTOR  Q 


EfTective  mobility  II  (Q-dependent  rip)  Components  of  the  mobility  (low  density) 


EQUIVALENT  OXIDE  THICKNESS  tgq  (nm)  EQUIVALENT  OXIDE  THICKNESS  tgq  (nm) 


Components  of  the  mobility  (high  density)  Mobility  vs  static  dielectric  constant 
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Experimental  confirmation? 

-  Poor  mobility  observed  almost  always,  but  mainly  due  to  poor  material 
quality  (traps,  interface  states,  non-stochiometry,  interfacial  layers). 
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using  Oguzman  et  aL  pseudopotential  parameters 
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Semiclassical  modeling  of  small  semiconductor 
devices 


Semiclassical  modeling  of  small 


Asen  Asenov,  Glasgow  University. 

Discrete  impurities  and  quantum  potentials  in  MOSFET  modeling 
Richard  Akis,  Arizona  State  University, 

Effective  potentials  for  quantum  effects  in  MOSFETs 

IBM  Research  Division 
T.  J.  Watson  Research  Center 
Yorktown  Heights,  NY  10598 

June  2001 
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Semiclassical  modeling 


Dopant  fluctuations 


•  Various  flavors: 

-  Strict:  may  use  the  bulk  (3D)  Boltzmann  Transport  Equation 
(BTE),  a  la  Kohn-Luttinger  (Liouville-von  Neumann  with  infinitely- 
may  randomly-distributed  impurities),  even  if  Quantum  Mechanics  is 
pervasive  (bands,  collisions,  etc.) 

-  Somewhat  looser:  may  use  a  lower-dimensional  BTE  (quantum 
confinement),  or  any  Master  equation  (dots).  No  transport  'In  the 
quantized  direction’ 

-  Very  loose  (‘psychologically’  classical):  Approximate  quantum  effects 
so  to  allow  the  concept  of  ’particle'  or  ’trajectory'  (Bohm  trajectories, 
quantum  potentials,  effective  potentials.  Wigner  paths/trajectories...) 

•  Anticipated  failures: 

-  Quantum  confinement  (strict):  channels,  wires,  dots 

-  Tunneling  (somewhat  looser):  band-to-band,  across  barrier,  resonant, 

-  Interference  effects  (somewhat,  looser):  Bohm-Aharonov,  current 
oscillations  in  FN-tunneling  and  channel  current  (?!),  ... 

-  2D/3D  confinement-effects  (somewhat-looser/very-loose?):  narrow 
channels  in  double-gate,  point  contacts,  strongly-coupled  dots 

-  High-energy,  high-field,  short-time  effects  (Somewhat /very  loose): 
Collisional  broadening,  intracollisonal-field,  coherent  phenomena,.. 

•  Virtues  of  semiclassical  modeling: 

-  We  have  survived  so  far... 

-  ’Intuitive’,  leading  to  deeper  understanding 

-  Simpler  (or  just.. .doable?),  allowing  more  accurate  description  of  effects 
which  maybe  dominant  in  m.iny  real-life  experiments/devices 

-  Determination  of  ’bulk  quantum  properties’  (band-structures,  strain, 
electron-phonon  coupling  constants.  Auger/impact  processes) 

-  Collective/many-body  effects  (Coulomb  effects,  quantum  or  MD) 


•  Important  (threshold  control  a  crucial  element  in  VLSI)... 

•  but  complicated: 

-  BTE  suspect 

-  Scatterers  closer  than  mean-free- path; 

«  no  rigorous  use  of  mobility  possible 

♦  multiple-scattering,  interference,  non  phase-breaking 

-  What’s  the  scattering  potential? 

•  So  far  limited  to  drift-diffusion  and  to  electrostatics  (probably  good  enough) 
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Quantum/efFective  potentials 


•  From  Wigner  or  Bohm 

•  ‘Phsychologically  easy’,  mathematically  hard 

•  So  far  approximated  (equilibrium)  or  of  a  heuristic  nature 

•  Of  high  appeal  because  of  their  'intuitive'  nature 


MVF  June  01  4 


o 

4^ 

•  pH 
00 
C 
cC 

C/D 

c 

o 

•pH 

5-h 

in 

+-> 

C 

U 

0) 

O 

a 

X-) 

c 

CC 

•pH 

TJ 

:3 

in 

C 


cC 

CD 


N 

c 


o 

■g 

> 

c 

D 

u 

CD 

CC3 

PQ 

u 

CD 

x: 

-M 

C 

CD 

D 

o 


rn  CJO 

.S^ 

C/D  cj 

C  CO 

o 

cj5 

C  i>. 

Si  ^ 

£ 

'O  c 

w  03 


CD 

4-J 

oo 

u 

•irH 

"c3 

'S 

E 

o 


"M 

CD 

c/) 

£ 

x: 

u 

< 


Si  ^ 

CO 

CJ 


3-^ 
5  « 

hi  00 

^-2 

s 

CD  CD 
-M  CO 


s? 

^  O 
CO  _rH 

CJ 

T3  o 
^  C 
,  CO 

^  -r-l 

CD  Xh  CO 

JS  >  cu 

■M  ^  U( 


ct 

2  o 

X  * 

CD  /]S 

C  o 

O  *rH 

.hn 

ti  ^ 

CO  N 
N 


D 

o 

x: 

CO 

Di 

C 

o  . 
T3  O 
CO  As, 
N  * 
Xh  H 
cq  u 
o  o 

O 

2  2 

^  s 

*0  J 


13  S 

8^ 


t-i 

6  CO 
S  DC 

'"'  DC  O 
CO  O  -S  ^ 

0)  J-I  o 

^  to  c  ^ 

KH  DC 


CO 

u 

-2  T3 

>  h3 

CO  CO 

CD 
X3 


T3 

C 

X 

£ 

CO 

CO 


tuo 

c 

•  rH 
^-) 

r2 

3 

CO 

C 


U 

,o 


CD 
£ 
CO 

o 

x: 

CO 

^  £  S 

•Th  CD  CD 

CO 

^  CO  CO 

QJ  CO  CO 

^  CU  o 

^  X  X 

CO 

•2  5  S 

-  °  I 

ui  «  ^ 

O  PQ  < 

<  I  I 


cO 

u 


DO 

c 


CO 

S  >-' 

Cu  Q} 

CJ 
CO 


C 

o 


D 

a 

(D 
U 
CD 

CO 

C 

o 

•  rH 
+-> 

CO 

C 

'a 

CD  C/J 

^  9 

g  ? 


C 

o 

-M 

CO 

3 

tj 

,2 

CO 

(D 

& 

3 

c 

CD 

£ 


CO 

C 

CO 

u 


u 

X 

a 


o  ^ 

-M  U 
cO  CD 
'X  N 


u 
CD 

U  't' 


c 

o 

6b 

C 

C 

CD 
CD 
u 
u 

a;^ 

CO  4_, 

C  u 

O  ,CD 


CD 

£ 


CD 


C3 

O 


(73 

§ 

Vh 

H 

;-4 

O 


o 

cri  P 

(73  (73 

O  pH 


cy^ 


cd 

(U 


(D 

H 


tS 

CO 


x:3 

a 

S) 

cd 

1 

K 

2 

•  pH 

XU 

m  pH 

g 

g  3^ 

4^  •i—t 

2  I 

<o  3> 
'O  ^ 

O  c 


N  - 

•S  ^ 


J 

t+H 


2 

<D 

.> 

s 

rv 

g 

PQ 

a> 

G 

D 

3 

o 

< 


Td 


o  c 
^  a 


CO 

CL) 

> 

•  pH 

5 


a  Q 


CO 

i 


I 

CO 

a 

•  pH 

I 

3 

CD 

£ 

CD 

4={ 

H 


cO 

a 

o 

N 


TO 

«-H 

•  pH 

PQ 

a 


a 

o 


CO 

D 

Dh 

o 


cO  c3 

00  .2 


CO 

*  ^ 
^  s 

4.^  CCS 

p£  T^ 

^  £ 
^  cO 

.£  CO 


I  i  •§ 
"  a 


CQ 


CD 

H 


I 


ClB 

5"  CO 
CJ 


TO 

C 

CD 


D 

CO  TO 

^  Si 

CO  C 

g-a 

•pH  r. 

+J  S 
CJ  ^ 
CO  . 

S  3 

DO  <-> 
C  t- 
O  0) 

h  C 
w  M 

U  > 

o  ^ 

P-r  CO 


c  ^ 

r/C 

co  2 

CD  C 

>  O 

>  .pH 

Q-h  cj 

°  2  o 

•a  -i 

E  C  H 

u 

2^  DO'S 
c  +j 
o  O 

s 

2  c 

8 

^  ^  -I 

.■y  2 

^  .-y  X 
C  c  o 
.£  -Cl 
O  Xh  CO 
U 

2  CD 

XP 

CD 

C 


CD  X 

CO  DO 

CD  -in 

x:  x: 

c 

CO 
CD 


T3 

CD 

E 

<4H 

c 

o 

CJ 


CO 

E 

<u 

•M 

CO 


o 

:§.2 


H  .2  53 


CO 

CO 

+J 

•M 

C 

c 

CD 

CD 

e 

E 

•pH 

•pH 

u 

u 

CD 

CD 

cx 

cx  . 

X 

X  CO 

CD 

CD  c 

o 

u 

•rH 

0) 

>,.2 

> 

CD 

U  a 

O 

X 

c 

o 


CO 

•d  « 

3 

o  u 

o*-i  o 

^  - 
^  C 


a 

£ 


X 

X3 

C 

o 

u 


CO 

a; 

£ 

CO 

u 

X 


CD 

N 


CD 
U 

C  S 

O  T3 

C  > 

<  w 

I  I 


The  metal-insulator  transition  in  open 
quantum  dots  and  arrays 

Jonathan  Bird,  Arizona  State  University  ^  ^  j  ^  • 

•  The  quasi-bound  nature  is  expected  to  give 

Much  recent  interest  in  the  observation  of  a  Kondo  rise  to  novel  signatures  in  transport  due  to 

effect  in  Coulomb  blockaded  quantum  dots.  electron  interactions. 
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Negative  electric  field 

respon-  sible  for  charge 

setback  from  the  interface. 
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the  result  is  not  symmetric  in  the  long  axis.  Thus, 
there  is  a  definite  input  side  and  a  definite  output 
side,  which  is  not  often  recognized. 


A  magnetic  field  mixes  the  states  and 

produces  an  entangled  state. 


Nanostructures  Research  Group 


Conventional  device  modelling:semi-classical 


Non-stationary  states  show  trajectories 


A  posteriori  quantum  hydrodynamics 


Quantum  potential  causes  repulsion  from  boundaries 
But  also  lowers  barrier  energies  to  allow  tunnelling 


Interpretation:  ,diSQQ¥er„atraje,Qtojry, 


No  overlap:  effective  statistical  mixture,  Overlap:  entangled  state, 

System  point  is  in  one  or  the  other  sub-state  Correlation 


2  fermions  in  a  harmonic  oscillator  potential 


FJow  near  singularities  The  quantisation  of  velocity 

circulation 


Traversal  Time  To  describe  vortex  motion  with  ab  initio 

„  quantum  hydrodynamics 

requires  Bohm  interpretation  _  f  v.dr  we  need  to  introduce  a  vector  quantum  potential  aQ(r,t) 


Compare  Lord  Kelvin’s  cigar  box! 
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Current  density  magnitude  transverse  across  the  cavity 


Reflecting  wall 


Incident  travelling  wave  in  high  transverse  mode 


Detail  of  upstream  density  and  velocity  flow  Upstream  flow:  full  form 
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A  posteriori  quantum  hydrodynamics 
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Application  to  MOS  device  modelling 
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Exponential  power 
spectrum 
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Coulomb  correlations  in  semiconductors  Coidomb  correlations  in 
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Dirk  Semkat/Michael  Bonitz  Semkat/Kwong/Kohler/Binder/MB 


Selfconsisicnt  evolution  and  decay  of  initial  correlations 


Dynamic  Structure  factor  of  interacting  electrons 

Solution  of  inhomogeneous  Kadanoff-Ba>aTi  Equations 


probe  photon  energy  (eV)  probe  photon  energy  <eV) 


Impact  Ionization  rates  in  semiconductors 


Mesoscopic  electron  clusters 

Easy  realization  of  strong  correlations 
Question:  possibility  of  Wigner  crystallization? 

Is  there  a  metal-insulator  transition? 


Radial  melting,  OM  -  Orientational  melting 


Single  particle  interference 

cease  to  effect  resistivity  rai  curve  ...  disappearance  of  the  weak  localization 
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EXPERIMENTAL  RESULTS 
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•  We  suggest  the  exponential  regime  indicates  a  transition  from  THERMALLY-BROADENED  to  •  While  our  simulations  show  features  reminiscent  of  the  exponential  regime  found  in  experiment 

ENERGETICALLY-DISCRETE  levels  in  the  dots  with  decreasing  temperature  they  are  UNABLE  to  reproduce  the  LOGARITHMIC  behavior  seen  at  much  lower  temperatures 


Phy*.  R«V.  B  fii  241302  (2001)  Phys.  R«v.  B  §2,  241 302  (2001) 
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•  IT  IS  OUR  HOPE  THAT  THIS  WORK  WILL  STIMULATE  RENEWED  THEORETICAL  INTEREST  IN 
THE  STUDY  OF  ELECTRON  INTERACTIONS  IN  OPEN  MESOSCOPIC  SYSTEMS 


-  B-field  induced  separation  of  +3/2  and  -3/2  hh  bands  in  k^space 

-  current  direction  determined  by  helicity  of  light  polarization 
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Elementary  definition  and 
properties 

Dynamical  equations 

-  Classical  force 

-  Infinite  potential  barrier 

Electron  -  phonon 
interaction 

Wigner  paths  and  MC 

simulation 

-  Paths  and  diagrams 

-  Multi  pi  Icity  of  Wigner  paths 

-  Quantum  self  scattering 

Two-time  Wigner  function 
f  (p,  co) 


Rossclla  Brunotti 
Paolo  Bordone 
Andrea  Bortoni 
Stcfano  Monastra 


-  Spectral  density 


THE  W][(GNER  FONCTION 


T he  Wigner-  function  approach  is  a  phase-space  formulation  of 
quantum  mechanics  that  allows  to  establish  many  analogies  with  the 
semiclassical  theory  based  on  the  concept  of  distribution  function 


Elementary  definition: 

P(riT2) 


Mr,P,t)  =  jdse-‘^'^^{r  +  s/2,t)'¥’{r-s/2,t) 


WIGNER  FUNCTION 


Properties  of  the  Wig-ner  function 


WIGNER  FUNCTION 

Potential  step  with  scattering  states: 


WIGNER  FUNCTION  FOR  A  DOUBLE  BARRIER 


COF3ERENT  PROPAGATIION  OF  THE  W.F. 


DYNAMICAL  EQUATION  OF  THE  W.F.  (i) 


// 


then 

Linear 

where 

/...('•.P)  =  \dse-“‘"'{r-\-sl2\<p,){<p^\r-sl2) 


From  the  Liouville  -  Von  Neumann 

with 


Hz:H„+V(r) 


Obtain 

V  ^  +  +  V(r),p]|r-s/2) 

at 


With  standard  elaboration: 

and 

Jrfse-''*'‘(r  +  s/2|X,[v(r),p]|r-s/2) 

=  \dp'V„{r,p-p')f„{r,p\t) 

where 

K{r^p)  =  ^^\dse*’"‘\y{r-sl2)-V[r  +  sl2)\ 


DYNAMICAL  EQUATION  OF  THE  W.F.  (2) 


Collecting  the  above: 


_9 

dt 


zA 

m  dr 


/».= 


For  "regular" potential  (up  to  quadratic): 


at  m  or  op 


as  in  Liouville  theorem. 


WIGNER  PATHS: 

The  iVigner  function  of  electrons  in  presence  of  potentials  up  to 
quadratic  evolves  as  an  ensemble  of  classical  particles:  each  point 
follows  a  classical  path  in  the  Wigner  phase  space  (Wigner  path) 


p 

Classical  phase  space 

r 

1 

A  1 

0 

1 

iL_ 

Quantum  mechanics 


P 

Wigner  function 

r 

Time  evolution  of  the  Wigner  Function 
of  a  wave  packet 


0  55  150  sec  255  300  3K>  «C0  500 


r[nni] 


Time  evolution  of  the  Wigne  Function  reflected  by  an 
infinite  potential  barrier 


WIGNER  P  A-THS  WIXM 
P  OX E NXI AL  S  C  AXXE  RING 

V  =  V+V'  F=-VF 

o  o 

iincar  or  quadraric  -  hnegra-dijferential  equation 

dt  m  or  op 

where 

F^(r,  p)  =  [f"  (r  +  s  /  2)  -  F'(r  -  s  /  2)] 

path  variahk's:  _ 

r(7)  =  r- p{t')lm  p{T)  =  p-'^r{t’)  dt' 

time  integraiior}: 

Lir,pj) 

as  in  Chamlx'rs  kinetic  integral  equation,  without  scattering  “out’',  with 
si  mi  far  inierprctalion: 

f(u 

— '4  m 

Neumann  expansion 


Wigner  ’’paths’': 

lUtJ  .  , , 

-  foj  rtf'i 

f  ’(0  +  . . . 

-  -4  ro) 


Boundary  conditions 


Phys.Rcv.B  58,  3503  (1998) 
Phys.Rev.B  59,  3060  (1999) 


Single  barrier  (4nm  x  .05eV) 


Wig-inier  IFiumction  confmedt 
by  .an  unifinilte  potential  barrier 


Obtained  by  means  of  W.  Paths  with  potential 
scattering  in  a  Monte  Carlo  approach 


\6 


For  the  infinite  bamerthe  Schroed.  eq.  holds  only  inside. 

For  a  box  (a.b) 

s(Jt)  _ 

r.ix,p,t)=  J  dse-"'^"‘^{x+s/2,t)'¥’{x-s/2,t) 


wliere 


!(,>;)  =  l-'(r  +  s/2) 


|2(/)-a-) 


if  x<{a  +  h)l2 
if  x>[a  +  b)ll 


17 


Wig-ner  Function  confined 
by  an  infinite  potential  lbarrier(ii| 


In  the  derivation  of  the  dynamical  equation,  the  derivatives 
at  tlie  boundaries  do  not  vanish,  and  we  have  an  extra 
integral  tenn 


+  i  ^dp'B  {x,  p-p')  £  y.  ( X,  p './ ) 


where 


B[x,p)=  sm 
m 


18 


SEPARATIION  BETWEEN  CLASSICAL 
FORCE  AND  Cyj[ANTUM  EFFECTS 


-,]dp'l\.{r,p-p')j\.{r,p\l) 


dt  m  dr 
where 

V,{r,p)  =  ^\dse“-"'p(r,s) 

( r,  5)  =  F  (r  +  5  /  2)  -  F  (r  -  5  /  2) 
or 

^A+P^  +  F^A  =  LJ^dp'VSr,p-p')fyr,p\t) 

dt  m  dr  dp  h  t 


where 


a^/, 

ar 


inJ 

j3(r,s)  =|^F(r  +  5/2)“^s/2 


f)F 

y{r~s/2)  +  ^sf2 

dr 


simulation  ... 


ELECTRON-PHONON  SCATTERING  WITH 
THE  WIGNER  FONCTION 

Extend  the  definition  of  Wipter  function: 

fy,ir,p,{nq}An\},t) 

=  2  +sl2,{nq }|p(f )|'- * ' 2,{«',  }) 

Hamiltonian: 

H  =  H,  +  +  \(r)  +  V'(r) 

where 

t2 

H  = - H  =y  b 

(r)  =  eEr 

V  fr)=  potential  profile 

Differentiate  the  definition  above  and  use 

if>j^P  =  l^,P] 

Obtain 

l/Jr.p.Kl.K,  M)  =  Z  X 

(r  +  i  /  2,  {n^ }  I  +  Hp*  +  +  V„(r)  +  V'(r)- p]|»- -  s  /  2,  {«', }) 


20 


ELECTRON-PHONON  SCATTERING  WITH 
THE  WIGNER  FUNCTION-  11 

Several  terms: 

^e-*-{r+s/2,{n,)\y,[H..pir-s/2.{n\])  =  -£'!f. 

M 

2e-'"(''+*/2.{«,)|>iK(r),plr-W2.(«',}>  =  -F^/. 

2«'*”('-  +  s/2.{«,)|>i[v'(r),pjr-W2,K,}) 

=  i^'V.(r,p-p') P' A”,), {»',}.!) 


Electron-phonon  term: 


Path  pariables\  time  integration,  Neumann  expansion  ... 


PHONON  AVERAGE 

Reduction  to  electron  Wigner  function 

Trace  over  phonon  variables: 

/„9p,tT=  X/wdp,n^SR,Si 
ncr 

Assume  separate  equilibrium  initial  density  matrix: 

/H.dp.R,Sn^S  =  oi~  f 

For  final  diagonal  termsy  with  initial  diagonal  W.f.,  each  Monte  Carlo 
term  contains  a  phonon  mode  twice  and  has  factors  like 

A/'wd  p,  s  i  ...cl,, + 1 1- ./wd  p,  R',  SR,  S)i 

The  simulation  must  be  thought  of  as  repeated  a  Uarge  number'  of  times 
for  each  {n^}.  Then,  if  hot-phonon  effects  are  ignored, 

for  absent  q: 

RJ 

^  for  real  or  virtual  absorptions: 

",  ", 

■=>  for  real  or  virtual  emissions: 


Plhioinion  scatteriinig’ 

Real  emission  ("in") 


Virtual  emission  ("out") 


Multiple  scattering 


''%hJ 

1 

r 

r — 

! 

IntracolUsional  fidd  effect,  interference  with  potential  scattering, ... 
MONTE  CARLO 
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Potential  step  -  one  phonon  scattering 


Multiplicity  of  Wig'jnier  paths 

Each  muhiple  path  can  be  obtained  with  different  diagrams 


i-r 

1^ 

Each  graph  contributes  with  a  factor 

Summing  up  all  terms,  the  contribution  becomes 

2cos((^.r.  -cot.j)2cos{q'r^ )  ->  2"]^  cos 


Double  barrier  -  one  phonon  scattering 


24 
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Qua  in  tu  m  self -sea  t  iteiri n 

IN  STRICT  FORMAL  ANALOGY  WITH  SELF 
SCATTERING  IN  SEMICLASSICAL  TRANSPORT 
SIMULATION 

Define: 

By  substitution  into  the  dynamical  equation,  an  extra  teim 
is  added  to  the  interaction,  proportional  to  F. 

Each  “free  flight”  in  the  Wigner  paths  is  multiplied  by  a 
factor 


corresponding  to  an  approximate  constant  life-time  of  the 
electron  state  (imaginary  part  of  the  self  energy). 

At  each  “scattering  evenf  ’  a  choice  is  made  between  self  or 
actual  scattering,  thus  correcting,  at  each  order,  the 
approximate  self  energy  with  the  exact  one. 


irWO^TIMlE 
WIGNEIR.  FUNCTION 


p(ri,ti;  r2,t2) 


=  ^'e'’®^^e“'^*^^(r+s/2|'Pd+T/2i)(4'd-T/2i|r-j/2) 
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TWO-TIME  WIGNER  F.  WITH  PHONONS  (i) 


H  =  H  +H  +H 

e  p  ep 


«K} 


Put 


Then 


TWO-TIMIE  WIGNER  F.  WITH  PHONONS  (2) 


Performe  traditional  transformation: 


t^L 

f  =  ^y^,  x  =  t,-t,  r  =  -^,  r  =r,-r^ 


r  +1* 


and  Fourier  transform: 

/..  (r,p,r,cd)  =  -ihG^  [r.p.Uco)  =  X  ] 

Kl  i 

Je"'‘'"i/Tg(r-?*72,{«J,r-r/2)g*(r  +  rV2,{/iJ,/  +  T/2) 
For  dynamical  equation  generalize  to  non  diagonal  phonon  states: 

/.('■>  p-  ®) = j  ' 

|e‘“’'2fTg(r-r72,{nJ,/-T/2)g*(r  +  r'/2,{n'J,/  +  T/2) 


and 


TWO-TIME  WIGNER  F.  WITH  PHONONS  (5) 


TWO-TIME  WIGNER  F.  WITH  PHONONS  (4) 


Dynamical  equation 

Je'"rfT|j^^g(r -r  72, {«,},/ -T/2)jg'(r  +  r'/2,(n',},<  +  T/2)  +  g^g'j 


If 


with 

=y  b/b 

p  La  ?  ?  » 

7 

then 

=  '■“({«,})«  +  ^(H.  +  H,  + 

Thus  the  frce-plKMion  dynamics  is  eliminated 


The  various  terms  in  the  dynamical  equations  are  treated  as  in  the 
standard  case: 

H  =-— V- 
2m 


yelds 


a? 


L{r,p,{n^],{n\],t,(o) 


V  =  F(r) 

velds 


a 


=  j  dp’  l\. {r.p-p')f,{r,p'.{n^},{ti'^},l.(o) 

where 


For  a  constant  or  harmonic  force  F 


d 

~dr 
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TWO-TIME  WIGNER  F.  WITH  PHONONS  (5) 


TWO-TIME  WIGNER  F.  WITH  PHONONS  (6) 


Taking  into  account  that 

h^.g[r-r'l2,{n^],t-Tl2) 

=  2,{„J|  b,.  I  'P  -T  /  2)) 

=  e""''("‘''yn,.+lg(r-r'/2,{/t,,...,n,  +  l,...},t-T/2) 
and  similar, 

the  electron-phonon  interaction 


yields 


I-'-®)  =  Z^(9'){ 


32 


The  general  equantion 

here 

I 


and 


A(gE'{''.}’{”',}-'-“)=E^(?’){ 

r»  *' 

7^/,(r,p  +  ft9y2, {/I,,..., -I, 

(r,  p  -  y  2, {/jJ \ '^.+ 1,..  j ,  ^ 

V^’''-“’'^^/^/^(r,p  +  ;i^y2,{/jJ, -l,...},^<y  +  tw,72) 


•  The  left  hand  side  of  the  above  equation  is  the  same  Liouvillian 
as  in  the  Boltzann  equation 

•  Path  variables  and  integral  form  of  the  transport  equation 

•  Neumann  expansion 

•  Wigner  Paths 

•  Monte  Carlo 

I"  At  eacli  interaction  “vertex”  half  phonon  momentum  and  half 
[  phonon  energy  is  trasfen  ed  to  or  taken  frotn  the  electron _ 
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Spectral  Density 
A(p.  <i))=fw-fp.6)yfw{  p) 


GaAs 

Polar  opcical  pt>onons 

T=300K 

up  to  3  scattering 

t=I50fs 


0  5  10  15  ao 


plAuJ 


Conclusions 


Wigner  paths 
he  defined  as  ballistic  | 
“flights”  separated  by 
scattering  processes 


Tbe  Wigrier-function 
approach  allows  a  description 
of  quantum  jmnsport  m  terms  of 
pihase-space  ^alo^s 
to  the  semiclassical  case  : 


A  Monte  Carlo  algorithm  can  be 
implemented  that  accounts  for  given 
initial  arid  boundary  conditions 


Momentum  and  frequency  can  be  considered  as 
independent  variables  in  a  two-thne  Green 
function  approach,  maintaining  the  path  Monte 
Carlo  approach 


PROB'-LEMS  - 


34 


35 


Density  Operator 

Time  Dependence/Liouville  Eq. 

A,,,(0=J]  I /(/)>  P(/) </(/)! 

I 

I  Coordinate  representation 

<  X  I  (/)|  x'  >  =  ^  <  X  I  /(/)  >  /’(/)</(/)  |x'  > 

i 


I  SRA 


ih- — ^ —  =  [//(/). /),,(/)  I 


f?<v| /){/)! \'> 

//; - 

(h 

f  tr  (  rT-  <r-  ^ 


Lm 


rJ\-  rl\'- 


■+r(x)-r(\')  <\|/?(/)ix'> 
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Density  &  Current  Density  Matrix 


p(x,x') ^<x\p„^,  (0|x'> 


Density  A(^)  X  \  (^)|  X  > 

jtx,x')  =  ^(V^-V^,)<x|p  (/)|x'> 
2nu _ 

Current  Density  j(x)  =  L//;?j(x,x') 
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Wnnost/ncti/ff'S  &  Tr.in^ivnK 


Real  Part  of  the  Density  Matrix,  T^O.O  K 
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IV 

Double 

Barrier 

RTD 


Current  vs  Applied  Potential  Energy  (77K) 
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RTD  as  a  Relaxation  Oscillator: 

(Verghese,  C.  D.  Parker,  E.R.  Brown  APL 
_ 1998) _ 


Transmission  Line 
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W.i/io.sffi/cf»r>'s  &  Trnn^ipnt^ 


Wc  will  concentrate  on  the  Wiancr  Function 


More  success  than  with  using  the  DM. 
Not  as  fast  as  the  DM. 


Wigner  Function/Equation  from  the  Density 
Matri.x  through  the  Wevl  Translbrmation 
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N.ini}‘itnictiiir>%  &  Tf-in^ionts. 
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The  Wigner  Equation 


Transient 
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Dissipation 


Diffusion 
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The  Wigner  Integral  Treated 
Analytically-Highlights  Correlations 
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Quantum  Mechanics;  The  Wigner  integral  represents  a 
correlation  of  slates  in  the  coordinate  representation. 
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V(x)-V(x’)  for  a 
single 

asymmetrically 
placed  barrier 


The  Five  Zone  Structure 


•  Five  zone  structure: 

--  The  end  regions  (two)  represent  the  contact 
regions; 

-  The  central  region  (one)--  incorporates  the 
quantum  barriers  and  wells; 

-  The  transition  region  (two) 

•  Variable  scattering  time,  smallest  in  contact 
region. 
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Device  Geometry 
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N.ino’itnicturf^  &  Tt,in^h'ni\ 


Structure  Emphasized 

The  Baseline  Study,  cant. 


C'onvcrged  tic  solutions  via  small  bias 
increments  show  no  hysteresis] 

[.arge  signal  transients  show  hysteresis. 
Relaxation  to  steady  state  oceurs  within  a 
pico-sec  and  is  bias  dependent 
Devices  sustain  steady  rcla.xation  oscillations! 
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Wigner  Function  NIN 


Wigner  Function  RTD 
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N.lnostnicrurt's  X  Trsnsinr 
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Current  Density  Versus  Applied  Bias 
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Applied  Bias  (ev)/kT 
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RTD  as  a  Relaxation  Oscillator 

(Lumped  Circuit  Elemenls,  RTD 
Represented  by  Winner  Equation) 


Wigner  Eciiialion 
Represents  the 
RID 


Realistic  r)c\  ice 
Response 
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d  Vp 
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dt 

Rr, 

di  _ 
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dt~ 

2;rZ„  1 

RO  Equations 

{'■(O-'/RC;)} 


Rr 


RO  Equations 
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Nnnn'iIrticUtrps  &  rf.ins/pfifs 
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RO  Equations 


di  _  Rp 


dt  IkT.,, 


^\il‘rUF.D  (0  -  *'o  (0  +  /  (0-^^  ’ 


v„  =  ("„/r,„  ;  =  ///,„  t  =  TlT„., 


k,  =  ^ZTcz,  =  27r JIc^,  R„  =  (■', / /,, 


SRA 
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RO  Equation  Parameters 


^0  ~  \/c  /  c  I, 

T,„  =27tjLC, 

R„  =  V,  /  /,, 


/  '  Observation:  For  a  fixed  Z,  oscillation 

characteristics  are  the  sanie  for  all  oscillatory 
periods  provided  the  N DR  is  not  dynamic. 
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RO  at  1 1.5  GHz:  Total  Current 
Density  and  Potential  Energy  Drop 

RO  at  1 1.5  GHz:  Particle  Current 
Density  and  Potential  Energy  Drop 
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Details  of  the  RO  at  1 1.5  GHz:  Particle 
Current  Density  and  Potential  Energy  Drop 


Details  of  the  RO  at  1 1.5  GHz:  Particle 
Current  Density  and  Potential  Energy  Drop 


First  part  of  the  oscillation 
Linear  I  versus  V.  L  R  ri^c  time 
Second  part  of  the  oscillation 

Drop  in  particle  current  and  increase  in  pulenlial  drop 
Third  part  of  the  oscillation 
-  Decrease  in  pt)iential  drop,  panicle  current  remains 
below  \ allies  assoeiatctl  with  tuM  part  tifthe  oscillation 
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N,ynn<;tnictiirp^  A  Tr-insirnta 


\  SRA 

N.ino’itnictiir 


‘Temporal  Invariance’ 

When  the  device  oscillates,  its 
characteristics  are  independent 
of  oscillatory  period! 


Potential  Energy  Drop  (nieV) 


I  SRA 

N.ino^truclutp’i  &  Transients 
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RO  at  1 1.5  GHz:  Total  Current 
Density  and  Potential  Energy  Drop 

RO  at  22.0  GHz:  Total  Current 
Density  and  Potential  Energy  Drop 
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N.inO'itiiictiirf'f  A  Tr.intiantK 


RO  at  41.8  GHz:  Total  Current 
Density  and  Potential  Energy  Drop 

RO  at  1 13.7  GHz:  Total  Current 
Density  and  Potential  Energy  Drop 
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Dynamic  ‘IV 


Potential  Energy  Drop  (tneV) 
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RO  at  1 1.5  GHz:  Dynamic  TV’ 

RO  at  22.0  GHz:  Dynamic  TV’ 

Potential  Energy  Drop  (nieV) 


Potential  Energy  Drop  (meV) 
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Total  and  Dynamic  TV’ 


Potential  Energy  Drop  (rneV) 
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Nnno^tnicliirp^  &  Tr,in<tipntii 


RO  at  6  CHz:  Total  and  Dynamic  TV’ 


RO  at  113.72  GHz:  Total  and  Dynamic 
TV’ 
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to^tritctiiros  A  Tran^.ipnt'i 


N.iiw^tnictiires  &  Tr,ir)f;iontii 


N.innstructiirt'^  &  Tr.in'iif’nt^ 


N,inoslnic1tirf''i  &  Tntnufpnts 


15  RO  at  22.0  GHz  (26  ps) 


N.irjostniCturPS  &  Tr,tns/rnts 
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NanostnscUiros  &  Transionts 
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Nano^tnictinfi’i  &  Tmnsipntu 


Nano^tniCtiirp'i  &  Trfmsipntfi 


N.inostructtirof:  &  Transients 


N.innstfticttires  &  Trans/rnts 


Nrinofifructurps  &  Ttrinsionts 


Potential  Energy  (meV) _ ^ _ Potential  Energy  (meV) 
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N.inn'itnictiirc’i  &  Transients 


Ntino’itructur/'s  &  Tr<ll 


I  RO  at  22.0  GHz  (50  ps) 


N.ino'^tntctiirp^  &  Tr.infiicnls 
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From  Simulations  (Post  Processed): 
Quantum  Potential,  Within  and  to  the  Left  of  the 
Barrier — For  Different  V^alues  of  Bias 


Conclusions  Suggested  by  the 
Quantum  Potential 


Simulations  are  consistent  with  the  filling  of  a 
well  defined  quantum  slate  prior  to  the  drop-back 
in  current. 

Quantum  potential  is  not  well  defined  after  the 
drop-back. 

Certainly  ihere  arc  coinpiilalitMUil  dirileultics  with  small 
\ allies  oTeliarge.  But  we  may  be  dealing  wilii  a 
mixture  of  a  large  number  ol\|uantum  states,  scattering 
Slates  and  dissipation. 
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Noise  Considerations 

•  For  nonlinear  systems  amplitude  noise  is 
not  an  issue  (M.  Lax  and  our  studies) 

•  Phase  noise  is  additive. 

•  Amount  of  phase  noise  depends  on  time  of 
disturbanee. 

•  Rccovei'y  from  disturbance  occurs  within 
one  cycle. 

•  Illustrate  with  ODE  and  Quantum  RO 


Injected  CniTent  Source  at  a 
Particular  Instant  of  Time 
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Same  type  of  current  source: 
different  time 


ODE/RO  Circuit  Response 
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ODE/RO  Circuit  Response 


Quantum  RO  at  1 13.7  GH/:  Fluctuation  Response 
of  Particle  Current  and  Potential  Energy 
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Conclusions 


Quantum  Van  der  Pol  oscillator  provides 
quantilalivc  values  ibr  the  large  signal  switching 
and  recovery  times  of  RTD. 

Large  signal  frequencies  near  120  GHz  for  non- 
opiimized  structure. 

Small  signal  oscillations  sustained  at  higher 
frequencies.  (f>670  GHz). 

Quantum  Van  der  Pol  oscillators  recover  from 
tluctuations  within  one  cycle. 

Phase  noise  dominates. 


Speculations 

Periodic  invariance  is  probably  due  to  time 
increments.  If  use  was  made  of 
femtosecond  increments  would  probably 
see  high  frequency  contributions. 

Quantum  potential  interpretation  not  as 
simple  for  higher  (>300  meV)  barriers. 


I  SRA 

M.mo'itfiicttirps  A  Tr.insirnl^ 


\  SRA 

Nnno'.tnictiirr'S  &  Tr.insir 


Potential  Energy  Drop  (meV) 
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Conventionai  TCAD 


Process-Device  simulation  (TCAD)  is  used  by  the  semiconductor  industry  to 
expedite  and  optimize  process  development 

TCAD  enables  the  efficient  use  of  experimental  cycles  and  provides  major 
cost/time  savings  during  early  design  iterations 

TCAD  device  models  are  based  on  semi-classical  transport  theory  which 
relies  on  approximations  to  facilitate  performance  efficiency 
Semiclasslcal  transport  models  will  become  increasingly  inaccurate  with 
decreasing  due  to  the  onset  of  quantum  effects 


0.25tim  0.1  Sum  0.1  Sum  0.1 3um  90nm  70nm  SOnm 

Experimental  calibration  will  be  increasingly  required  to  calibrate  TCAD 
device  simulators  which  will  undermine  their  usefulness 
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Comparison  of  Contemporary  Device  Simulators 


Nonequilibrium  Quantum  Field  Theor 


Poisson  Equation 


Dyson  Equation 

[r-w-i’jj  = ; 

Kinetic  Equation  | 


Current  (l-V) 
CapacitancB  (C-VO 


NQFT  is  the  most  comprehensive  theoretical 

framework  for  device  simulation 

Within  a  single  numerical  platform,  NQFT  can 

encompass 

-  Quantum  transport 

-  Scattering 

-  Atomic-scale  material  properties 

The  semi-classical  Boltzmann  equation  is  an 
approximate  form  of  the  NQFT  kinetic  equation 

Boltzmann  Equation 

|-/(v,r,0  +  vV/(v,r,/)  +  — FV5/(v,r,0 
dt  m 

=  Xl/(v'.r.05(v.v)-/(v,r,05(P,v)] 

NQFT  Kinetic  Equation 

“  Co/nowfat/on»/ 


Principal  Equations 

[x  -  )  - i' (£,  /fcj. )]  G' (£, =  /  ■“ 

Ctj(.EA)  =  Sc,' + 

LI 

Scattering 

t“(£)  h'-  ‘ 

Scattering  is  both  the  key  and  impediment  to 
numerical  quantum  transport  simulation 


-  Impulse  response 

-  Convolution  integral 
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•  Definitions 

H^(kJ  -  Tight-binding  Hamiltonian  goC^'.^x)  -  Zero-order  correiation  function 

t  *  (£,  )  -  Retarded  self-energy  (out-scattering)  Self-Energy  (in-scattenng) 

G *(£,)[  )  -  Retarded  Green  Function  (Propagator)  G^j{E,k_)  -  Full  correlation  function 

Kj  (£,  k,,E\k[)  -  Interaction  Potential  (Green  Function)  hO^ 

'  Computational 

-  Physical  Sciancas  Research  Laboratory  ■■  Materials  Group 


Noneauilibrium  Quantum  Field  Theory 
Physical  Observables 

Density  of  States 

{quantum  mechanical  DOS) 

Electrons 

=  //(£, 0^)4,  (£,  it  J  (equilibrium) 

=  i/(£,  (£,  )  +  f/(£'.^*)/4*(£,iFx) _ (non-equilibrium) 


fdEr  dk,  ^ 

n.  =  -2i  I  —  1  — (£,  k^ ) 

'  J  2;r-' (2;r)' 


(electron  density) 


^,  =  J  J  ^  [w,.,.,  (*x)G, ‘x)  -  C,‘ .1  (E.  tx  w.,x(*x )]  (electron  current) 

Holes 

G,’(£jj^)  =  -i[\-  fiE,(f>p)]A..{E.kJ  (equilibrium) 

=  -i[l-/(£,0^)],4,^(£,it^)-i[l  -/(£,0*)}4'(£.jt^)  (non-equilibrium) 


,frf£f  dk,  - 

Pi  -  “2' J  ^  (hole  density) 

°  f  1  [w„..  (*X  )C;.,x(E,  *x)  -  „  (E,  *x)W,.u(ix )]  (hole  current) 


J  2;r-'  (2;r)" 
2e  rdEf  dk,  r 

h  J  2;rJ  (2«-)' 
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Nonequilibrium  Quantum  Field  Theor 

Spatial  Discretization 
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Block  Tri-Diagonal  Hamiltonian 
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it  e  l  t 
ltd  t 

i  t  e  ...  t  j  . 

:t  e  t  t 
l  let  t 
t  t  e  t  t 

.  t  .  t.  £ . t  . 

t  e  t  :t 

t  t  E  t  M 

t  t  e  t  I 

. t . t.  £  ^  . t  . 

t  e  t 
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Presently  using  multi-valley  effective-mass  based  discretization 
Simulated  bandstructure  is  intrinsically  linked  to  grid  morphology  in  lateral  (z) 
direction  ^ 

E(k^)  =  -^ — [l“COs(A:.Aj] 
w.A, 

Uniform  grid  is  required  for  the  lateral  (z)  direction  to  avoid  spurious  reflections 
Non-uniform  grid  can  be  used  in  vertical  (y)  direction 
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Nonequilibrium  Quantum  Field  Theory 

Boundary  Conditions  for  Green's  and  Correllation  Functions 

Boundary  condition  assumptions 

-  Contacts  are  modeled  as  translationaliy  invariant  open  boundaries  and  are 
assumed  to  be  in  local  equilibrium 

-  Closed  boundaries  are  set  where  and  G*'  can  be  safely  assumed  to  be 
negligible 

•  Boundary  eigenvalue  problem 

[z.+{zj"-^k.=0 


Boundary  Green’s  function 

K^=QZQ''go\ 

SL  =[E-i-iezQl' 
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Nonequilibrium  Quantum  Field  Theory 

Recursive  Solution  of  Dyson  and  Kinetic  Equations 

General  Dyson  equations  for  O'*  and  G* 

G,y=g,' 

J.k  J.k  l.k 

Recursive  Greens  Function  (RGF)  Technique 

-  Direct  solution  of  the  Green’s  functions  is  computationally  prohibited  (NjNy)^ 

-  Physical  observables  are  clustered  along  the  Green’s  function  diagonals 

-  Dyson  equations  can  be  used  to  recursively  construct  Green  and  correlation 
function  diagonals  in  a  two-step  manner 

-  Perturbative  elements  are  the  off-diagonal  elements  of  the  Hamiltonian 

y.* 

st/  =  Sh  +  S  +  E  gfjhkSlf 

M  J.k 

-  Self-consistent  scattering  can  be  efficiently  Incorporated  into  the  RGF 
approach 
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Nonequilibrium  Quantum  Field  Theory 

Recursive  Solution  of  Semi-infinite  Green/Correlation  Functions 


Nonequilibrium  Quantum  Field  Theory 

Recursive  Solution  of  Full  Green/Correlation  Functions 


Sii 

9  o  a 


si 


E3>  .  . 


SiJ  8i,i  8i.i^i.i-t\8i+],i 

~  8>.i  8iJij*[8Hu*\^i*Ui8ij  8iJij*[8i^ij+\^i+\j8ij  8iJi.M8i*\jHh*\j8ij 

8/!=[(s.’r-ut«'u,tul' 

8ij  ~ 


Physical  Sdancas  Rasaarch  Laboratory 


Computational 
■  Materials  Group  CJ^ 


/^/l  _i_  7T ^ 

Gij^8u  +8ui 

Uiiu 

j+is»+i,i+r 

J 

Ou^gu  +gu^i. 

uul"i-^8i.ni 

Hi.i8u 

Physical  Sdancas  Rasaarch  Laboratory 


Computational 


Noneauiiibrium  Quantum  Field  Theory 

Parallel  Implementation  of  the  Recursive  Green’s  Function  Method 

Parallel  energy  integration 

-  Green’s  function  calculations  at  each  energy  are  distributed  to  slave  processors 

-  Found  to  be  highly  inefficient  for  scattering  simulations  due  to  varying  execution 
time  per  node 

Parallel  Recursive  Green’s  function  algorithm 

-  Spatial  grid  topology  is  distributed  across  processors 

-  Cross  element  coupling  terms  are  calculated  in  parallel 

-  Boundary  Green’s  functions  are  calculated  and  disb'ibuted  to  enable  parallel 
calculation  of  G*^ 

-  Correlation  functions  are  then  Iteratively  calculated  until  scattering  self- 
consistency  Is  achieved 


-  80%  parallel  efficiency  is  exhibited  for  60  processor  simulations 

ComputationMl 
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Scattering  Model 


Kinetic  Equation 


Current  Continuity 


|rfr,[r(i,jc,)C*(x„x')-r(x,Jc,M(x„x)]=0 

a 

Rate-Based  Scattering  Model 

Approximation: 


-  Enables  efficient  self-consistent  scattering  calculation 

-  Simple  calibration 

-  Accoi/nte  for  resistive  effects  (e.g.  L^ff  scaling) 

-  Physical  SciatKas  Raaaarch  Laboratory  . -  U 


Computatiorral 
Uatarials  Group  CJ^ 


Quantum  Mechanical  Simulation  of  MOSFETs 
Numerical  Issues  and  Computational  Platforms 


•  Algorithms 

-  Numerical  NQFT  solutions  rely 
heavily  on  matrix  inversion  and 
multiplication 

•  Computing  Platforms 

-  Commodity  clusters  are  the 
enabling  platforms  for  2D  NQFT 
simulation 

-  All  runs  were  performed  on  a  110 
node  2xPII!  450MHz  Linux  cluster 
which  cost  under  $250,000 

>  Performance 

-  90nm  Ballistic  simulations: 

-1  hour/bias  point 

-  90nm  Scattering  simulations: 

-8  hour/bias  point 
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Waveguide  Resistor  Transmissiort  Probability 

. .  . . I  ■!  .  t  ,  ,  , 

B- - Ballistic  p;^- 

- L=100nm  <2!:=* 

,  ‘  V=0.5V  r 


Empirical  Scattering  Model 

Resistive  electron  wave-guide 


Energy  (eV) 

Waveguide  Resistor  Transmission  Probability 


'  - L»1nm 

■ - L-IOnm 

1  - L=100nm 

:  V=0,5V 

■  1  I  1  1.  1  .l.l 

0.4  0.5  *  Ballistic  Case 

-  Subband  structure  is  apparent  and 

m  Probability  retained  up  to  1 0Onm 

-r-r-.-r-r-r^  -  Variation  in  L  causes  subband  smearing 

J  but  no  change  in  transmission  envelope 

-  Resistive  effects  are  not  accounted  for 

>  •  Constant-Rate  model 

'  Accounts  for  incoherence  effects  and 

back-scattering 

-  Subband  structure  disappears  for 
L=100nm 

-  Attenuation  of  transmission  envelope 
_ occurs  for  increasing  L 

-  Can  accounf  for  resisiiwe  effecfs  with 

‘  proper  calibration  ^ 

Computational 

Physical  Sclancas  Raaaarch  Laboratory  '  Uatarials  Group 


-|  MOSFET  Quantum  Mechanical  Effects 
Sub'bands 

•  Intersubband  scattering  processes  facilitate 
carrier  capture  under  drain  bias 

•  Back-scattering  is  effectively  suppressed  In 
the  post  source  barrier  region  by  coupling  into 
high-barrier  subbands 

•  As  a  result,  scattering  In  the  vicinity  of  the 
source  barrier  dominates  transport 

•  This  effect  is  captured  in  our  empirical 
scattering  model  despite  the  lack  of  relaxation 


MIT  90nm  Scattering  vs  Distance 

■  T  T  f  1  'f  «  I  1  I  I  '  I  I  »  I  I  I  I  ■' 

;  - Iniertace Bulk  Scattering 

L- - Bulk  Scattering 
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Spatial  Extent  of  Scattering  (Ralativa  S/0  Diatence 
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Quantum  Mechanical  Effects  in  MOSFETs 
Subbands 

MU  SOnm:  Tranamlaaton  Probability  _ MU  90nm:  Current  Density 


-r=1moV  (1.52x10*V) 

"  r=5m«V  (7.60x10”f ') 
••r=12.5mBV  (1.00x1  o’ V)  j 


“r-lmaV  (1.52x10’ V) 

"  r=5m«V  (7.60x10"s’) 
“r»12,5meV  (I.OOxIO’V 


0.2  -0.1  0  0.1  0.2  0.3  0.4  O.S  0.6 

Energy  (#V) 


0.1  0.2  0.3  0.4  O.S 

Energy  (eV) 


Ballistic  transmission  probability  shows  steps  at  the  onset  of  primed  (4-fold 
degenerate)  and  unprimed  (2-foId  degenerate)  sub-band  conduction 
Spatial  localization  of  T  induces  coupling  between  sub-bands  which  smears 
the  transmission  characteristics 

Increase  in  r  causes  decline  in  transmission  envelope  due  to  back-scattering 


Phyvcal  Sciancas  Raaaarch  Laboratory 
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Quantum  Mechanical  Effects  in  MOSFETs 
Charge  Centroid  Shift 

Charge  centroid  shift  Is  the  predominant  quantum  effect  in  MOSFETs  with 
highly  scaled  t^^ 

Inversion  layer  charge  shift  increases  effective  t^^.  reduces  induces  V,  shift 
and  degrades  sub-threshold  performance 

Source/Drain  centroid  shift  induces  the  formation  of  an  interface  dipole  with  a 
couesponding  charge  build-up  to  satisfy  macroscopic  charge  neutrality 


Channel  Electron  Concentration 


Source  Concentration 


0  2  4  6  8  10  12  14 

OlsUnca  (nm) 
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Quantum  Mechanicai  Effects  in  MOSFETs 
Charge  Centroid  Shift 

Source/Drain  centroid  shift  induced  dipole  creates  a  large  electric  field  along 
the  extension  interface 

Cuuent  density  is  tightly  distributed  along  the  interface  implying  that  cuuent 
crowding  may  be  exaggerated  in  conventional  series  resistance  calculations 


ional 


Physical  Sciancas  Rasaarch  Laboratory 


Computational 


Quantum  Mechanical  Effects  in  MOSFETs 
Ballistic  Sub-band  Mixing 

>  Schrodinger  Equation  in  Spatial  Representation 

(v;  +  v; )+  £,(.v,  z)]r  (v,  z)  =  EV'C.V.  z) 

*  Discretized  Schrodinger  Equation  in  Mode  Representation 
(£-£,  -£’>?-  =  0 
rr  =^idyt>r(y)r,(y) 


Double-Gate  MOSFET 


Bulk  MOSFET 


Modes  are  de-coupled 


Modes  are  coherently  coupled 
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Intrinsic  Mode  Coupling  in  MOSFETs 


Double-Gate  10nm/2.5nm  L  MOSFET 

.5  - - - ,  ,  ,  , - I  ,  ,  ■■  ■■  i  0.5 


90nm  L  Bulk  MOSFET 


4  S  6  7  8  9  10 

Subband 


1  2  3  4  5  6  7  B  9  10  11  12  13  14  15 
Subband 


Double-Gate  MOSFETs 

•  Minimal  mode  coupling  occurs  due  to  uniformity  of  channel 

•  Well  approximated  by  de-coupled  subband  models 
Bulk  MOSFETs 

•  Significant  intrinsic  mode  coupling  into  high-order  modes  occurs 
under  on-state  biasing 

•  Effect  is  not  captured  by  semiclassical  transport  models 
incorporating  sub-bands 

-  Physical  Sciancas  Rasaarch  Laboratory  . Uatarials 
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Quantum  Mechanical  Effects  in  MOSFETs 
Sub-Threshold  Tunneling 

Sub-threshold  tunneling  through  the  source  barrier  enhances  off-state  leakage 

and  compromises  ttie  lon^U  trade-off 

50nm  MOSFET  exhibits  a  -25%  l^ff  Increase 

25nm  MOSFET  exhibits  a  -40%  increase 

Further  complications  may  arise  from  trap-assisted  tunneling  due  to  S/D 

extension  lateral  dopant  diffusion  into  the  channel 


Normalized  Transmiasion 

t  (  ■*  »  r  I'l  I  I  I  I  I  r  ■  I  I  I  ■  I  I  I 
- MIT  Mnm 

- MIT  Mmn  Jj 

- MITJSnm  (V  -1.010  /i 


Normalized  Current  Density 
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Si  Spectral  DOS 

Phosphorus -doped  vs  urxJoped 

— ' - r~~ - ' - 1 - - - 1  ' 

Intrinsic  Si 

P  doped:  £o=45meV,  eV*.T=6.6x10’’®sec 

P doped:  £p=45meV,  eV*,t=6.6x10"’*sec 


-1  -0.5 

Energy  (eV) 
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Quantum  Mechanical  Effects  in  iVIOSFETs 
Id'^gs  Saturation 


Double-Gate  SOI  Simulation  Study 
Saturation  of  Gate  Control 


l-G  Simulation  Boundary 

Double-Gate  MOSFETs  are  predicted  to  display  optimal  scaling  properties 
and  carry  the  potential  of  realizing  (sub?)  10nm  technologies 
However,  electrostatics  dictates  that  LQ/Tg,  remain  roughly  constant  to 
suppress  short-channel  effects  and  achieve  target  ‘on/off 


•  Set  Le/Tgj  =  4  to  control  SCE  ...  keeping  tox=1.5nm  constant 
>  Evaluate  device  performance  for  10nm  <  Lq  <  80nm 
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DGSOI:  lOnm  L„,  2.5nm  T., 


~NQFT 

-  DESSIS  (Drift^DiffusionJ/ 


QUARTZ  Simulation:  T.,=6nm 


Loss  of  " 

Source  Barrier 


0.2  -  //  Unstable  DESSIS:  « 

yV  self-consistency  ■ 

0 , :  ■ , .  1 .  . .  ■ .  1 . .  ■  1 . .  .1  .1.5 1 

0  oV  04  0.6  0.8  1  i;2  "n..  *  " 

y  Distanc*  <nm) 

'  Drift-diffusion  simulations  show  an  instability  in  self-consistency  for  Ves>0.8V  due  to 
breakdown  of  macroscopic  charge  neutrality  at  the  source  boundary 
'  Boundary  charge  neutrality  is  maintained  for  quantum  transport  simulations  by  Invoking 
the  equilbrium  approximation  in  the  S/D  extensions 
'  Quantum  transport  simulation  predicts  current  saturation  with  Vcs=0.8V  due  to 
elimination  of  source  barrier 
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Double-Gate  SOI  Simulation  Study 
Impact  of  Current  Saturation  on  Device  Scaling 


NQFT: 


Drift  Diffusion:  L  fT  »  4 

■  I  ■  ■  ■  I  .Ti  I  P  I  y  I  rr 

- L  ■lOnm,  T  »2.5nm 

a  ’  II 

—  Lg=20nm,  Tj=5nni 

. L  »40nm,  T  >10nm 

■  -  L  =80nm,  T  =20nin 


0.4  0.6 


0.4  0.6 


'  Cases  studied  have  LgH's,  =  4  to  maintain  control  of  SCE  with  Lq  scaling 
-  QUARTZ  predicts  l^gg,  @  Vqq^I.OV  increases  with  decreasing  Lq  down  to  ~30-40nm 

•  Subsequent  reductions  in  Lq  lead  to  reduced  l^^g,  @  Vqq=1  .OV  due  to  elimination  of 
source  barrier  and  current  bottleneck  with  decreasing  T^ 

•  A  technology  with  ^OD  =0.6V  show  increasing  l^^g,  with  Lq  for  all  cases  studied 
>  DESSIS  shows  reduction  of  Ipsai  at  LQ=10-20nm 

Computational 

-  Physical  Sciancas  Rasaareh  Laboratory  . Materials  Group 


►I  Transmission  Characteristics  DOUbiG-GatC  SOI  SimUlatlOn  StudV 


•  Decreasing  T,j  is  predicted  to  enhance 
mobility  by  breaking  valley  degeneracy  and 
preferentially  occupying  low  m.^  valleys 

'  Macroscopic  charge  neutrality  must  be 
satisfied  in  the  source 

•  Velocity  enhancement  is  cancelled  out  by 
decrease  in  density  of  states  in 
transmission  characteristics 

Decreasing  T„  degrades  max 
...  cannot  scale  indefinitely 

Eiectron  Oansity  vs.  Body  Thickness 


1  0  0.1  0.2  0.3  0.4  0.5 

en«By  (eV) 

DGSOI  Current  Density 
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Quantum  Transport  Simulations  of  MOSFETs 
Scattering  Model  Calibration 

Mrr  90nm:  Quantum  Ballittic  v»  Experimental  I  MTT  90nm:  Quantum  Ballistic  vi.  Experimental  I 

10"*  f  »  »  »  I . .  . . .  0-35  I  . .  I  1  I  I  m  I  .  I  1  . I 


-  Experimental 
-Quantum  Ballistic  Limit 


«■*  r  v„«i.oiv 


-  Experimental 

-  Quantum  Ballistic  Limit 


MIT  90nm  Device 

Quantum  Mechanical  vs  Semiclassical  Transport 

MIT  90nm:  Quantum  vs  Semiclassical  MIT  90nm:  Quantum  vs.  Semiclassi 

*  »  . .  '  r  0.1  ITT,  M  »  »  »  »  1  '  ’  '  »  I'TTl-l  fl  TTTfTT 

V  “1.01V  ■  *  Experimental 

‘  r  “  ’  ^  ■  - Quantum  (r=12.5meV) 

jFy''  0.00 - Drift-DIffusion 

•  .  '  \  Hydrodynamic 
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Discrepancy  in  quantum  mechanical  sub-threshold  current  is  due  to  non-locality  of 
self-consistent  charge 

Device  drive  current  is  25%  of  theoretical  limit 
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Experimental 

—  Quantum  {r=12.5meV) 

—  Drift-Diffusion 
■•••  Hydrodynamic 
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Hydrodynamic  model  shows  minor  current  discrepancy  in  sub-threshold  regime 
Drift-diffusion  model  shows  minor  current  discrepancy  for  saturation  biases 
Quantum-mechanical  simulation  using  single  parameter  rate-based  scattering 
model  shows  good  overall  fit  with  experimental  data 
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CMOS  Effective  Mobility 
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EFFECTIVE  FIELD  [  MV/cm  ] 


effective  FIELD  C  MV/cm  ] 


Tak^.  at  ai.  IEEE  Trana  Bectron  Oavicae  (41)  12.  2257  {1994} 

>  Effective  mobility  is  a  measure  of  the  local  mobility  averaged  over  the  vertical 
carrier  density  in  the  channel 

'  Both  nMOS  and  pMOS  show  a  decrease  In  effective  mobility  with  increasing 
effective  field  (inversion  charge  density) 

'  Roll-off  of  mobility  has  been  generally  attributed  to  interface  roughness 
'  Atomistic  interpretation  has  been  lacking 
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Non-Equilibrium  Quantum  Field  Theory  Calibration 

Mobility  Model 

•  Linear  Response  Conductivity 

^*^1  ~~  I  (Current-Current  Response  Function) 

•  Conductivity;  Final  Form 

^  h  ^  dE  ^  2n  I  v„r„(£)  vl 

•  Effective  Mobility  and  Corresponding  MOSFET  Test  Structure 

< -  LQ«-»oopm  - ► 


Y‘  ^ 


Vos-OmV  • 

Vqs  modifies  Computational 

Physical  Sciancas  Rasaarch  Laboratory  ■  Maferfafs  Group 
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Non-Eauilibrium  Quantum  Field  Theory  Mobility  Model 

Single  Parameter  Scattering 


NQFT-1D  Simulation  of  Takagl  Device 


QUARTZ  Simulated  p  vs  Experimental 
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•  QUARTZ  M^»7,2x10'‘cm’ 

ZL 

■  QUARTZ  N^-3.0»10"em’’ 

SIO,  .  , 

SI 
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♦  QUARTZ  N^»7.7x10''cm'’  \ 

Distance  (nm)  (MV/cm) 

Single-parameter  scattering  model  brings  mobility  into  empirical  range  for  r=8meV 
Univsersal  mobilitv  behavior  is  exhibited  with  the  NQFT  mobility  model 
Strong  disag rement  with  empirical  universal  mobility  curve  begins  to  occur  for  high 
effective  fields  >  0.3  MV/cm) 

Single-parameter  scattering  model  alone  is  insufficient  for  capturing  effective 
mobility  behavior 
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Quantum  Transport  Simulated  Mobilitv: 
Bulk  +  Interface  Scattering  Model 


QUARTZ  Simulated  p  vs  Experimental 


- S, 

f-Co  +  iy 


Parameters:  ^ 

,  I 

A„  =  I.SeV  (coupling  strength)  .! 

e  D  =  0,42  (defect  energy  above  E J  = 

y  =  -  =  90meV  linewidth  (lifetime) 


"  QUARTZ  N^0.0«10' 

*  QUARTZ  N  "T.Tjiia' 

*  Qif;yTTj:« 


Interface  mode!  assumes  a  strong  peak  of  localized  states  exists  in  the 
Si  conduction  band  due  to  atoms  in  the  Si-SiO^  transition  layer 
Qualitative  view  of  the  interface  suggests  that  there  should  be  a 
transition  layer  that  averages  Si  and  Si02  conduction/valence  band 
properties 

Bulk  +  interface  scattering  model  produces  universal  mobility  behavior 
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Si-SiO,  Mobility  Model  Validation 


■  Direct  experimental  validation  Is  difficult  due  to  mixing  of 
interface  and  bulk  states 

>  Density  Functional  Theory  (OFT)  electronic  structure 
calculations  can  provide  Information  on  structure  and 
density  of  states 

Existing  DFT  work  has  not  focused  on  the  detailed 
interface  properties 


DFT  Calculations  (Lfu  and  Stumpfi 

•  Performed  using  VASP  on  CMG  Linux  cluster 

•  Examined  three  interface  structures: 

-  Si  -  SiO]  (P^quartz) 

-  Si  -  SiOj  (trydimite) 

-  Si  -  SiOj  (^rystalobalite) 

*  Two  basic  interface  Si-O-Si  bonding  structures  observed: 

.  Si'O-Si  Bridging  bonds 
.  Si-O-Si  Extended  bonds 

*  r/je  bridging  and  extended  interface  bonds  are  present  in  all 
5/O2  structures  and  are  periodically  distributed  along  the 


SI-O-SI  ^(*1 
bond 


VS. 
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CalibratBd  localized  interface  scattering  model  shows  poor  agreement  with 
both  experimental  data  and  drift-diffusion  simulation  for  high  Vq^ 


(miW) 


Interface  scattering  may  depend  on 

-  Non-tetrahedral  atomic  arrangement  of  bridging  and  extended  interfacial  0  atoms 

-  Distortion  of  Si  atomic  positions  below  the  Si-SiOj  interface 

If  interfacial  O  and/or  distorted  Si  atomic  positions  are  responsible  for  high-field  {!.„ 
behavior,  there  should  be  some  correlation  to  Si  areai  density  (p^j) 

Takagi’s  data  shows  that  (110)  Si  orientation  should  have  the  strongest  interface 
scattering  followed  by  (111)  then  (100) 


Takagi.  Mai.  IEEE  Trana  Electren  Omicaa  (41)  12. 2363  (1994) 

2  o  py  4 

p„(100)  =  ^  p(ll0)  =  ^  Ps,(ll0)  =  -^ 

o  a  \jo 

=  1.41p^.(100)  =1.15p^,a00) 

Scattering  from  tnterfacial  O  andlor  distorted  Si  lattice  as  the  cause  of  interfacial 
degradation  is  supported  by  experimental  data  for  orientation  dependence 

Computathnal 
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Atomic  Interface  Scattering 


'  As  Vq5  increases,  carriers  are  increasingly  driven  into  the  interfacial  0  layers 
where  there  is  a  modification  to  the  crystalline  periodicity  and  DOS  leading  to 
coherent  back-scattering 

>  As  Vqs  increases,  carriers  are  also  increasingly  driven  into  the  region  with 
interfacial  Si  distortion  where  there  is  additional  modification  to  the  crystalline 
periodicity  and  DOS  leading  to  coherent  back-scattering 

>  Long  range  distortion  of  the  interfacial  atomic  positions  corresponds  to  interface 
roughness  which  may  also  play  role ... 

ncoherent  Picture 

>  Phonon-assisted  trapping  in  the  interfacial  O  atoms  de-phases  electrons  and 
leads  to  Incoherent  backscattering 

>  Other  mechanisms  ?? 


>  Model  is  supported  by  Si  orientation  studies  on  behavior 

>  Anomalous  enhancement  of  nMOS  mobility  in  strained  Si  devices  may  be  explained 
by  modification  of  interfacial  atomic  position  leading  to  enhancement 
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Quantum  Transport  Simulated  Mobility: 
Atomistic  Interface  Scattering  Model  Calibration 


Quantum  Transport  Simulated  Mobility: 
Atomistic  interface  Scattering  Model  Calibration 


Interface  scattering  model  calibration: 

•  A  monotonicaily  increasing  is  used  to  account  for  the  increasing  interfacial 
O  and  distorted  Si  DOS 

•  The  Interface  model  is  spatially  extended  O.Snm  into  both  the  Si02  and  Si  to 
account  for  the  spatial  extent  of  the  interfacial  0  and  distorted  Si 

•  Model  is  calibrated  to  effective  mobility  data 

NQFT  Slmulatod)!  vs  Experlmont 


‘  Incorporation  of  the  effective  mobility  calibrated  interface  model  leads  to  excellent 
theory-experiment  agreement  for  2D  MOSFET  quantum  transport  simulation 
>  Quantum  transport  MOSFET  simulation  can  now  serve  as  the  basis  for  TCAD 
device  simulator  calibration  and  device  engineering  studies 


MiT  90nm  MOSFET  I 


MIT  90nm  MOSFET  I  -V^ 

. 

i  ExpsrtmMit 
—  Quantum  Transport 
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~  TakagI  300K  ^ 
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Exparlmant  1 

-  Quantum  Transport  -j 

-  Drift-DiffusJon  > 
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Quantum  Transport  Simulated  Mobility: 
Ballistic  Transport  vs  Scattering 


>  Ballistic  quantum  transport  simulations  reveal  that  the  MIT  25nm  device  is 
operating  at  -30%  of  the  ballistic  limit 

>  Interface  scattering  dominates  25nm  MIT  device  Iq,^ 

>  If  the  atomic  interface  scattering  model  is  to  be  believed,  there  is  little  hope  of  ever 
achieving  ballistic  performance  in  MOSFETs 


2Snm  Quantum  Transport  Simulations  MIT  25nm  Quantum  Transport  Simulations 


MIT  MOSFET  Simulation  Comparison 
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MIT  50nm  MOSFET  I 
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MIT25nm  MOSFET  I  -V^, 


mTr*niport(V^,M,CV)  j 

-j 
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MIT  50nm  MOSFET  t-V„ 

—  Qtanium  Transport  (V^ol.ZV) 

—  Qutmum  Tr»n»port 

- DfW-DtXuUon  (V^  JV) 

— •  Ot«.&fludon  (v'  -0  tV) 
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MIT25nm  MOSFET!  -V,. 


—  Quontuip  Tran*p»rt  (V^»1< 
CluinluRi  Tranfport  (V^MI.1 

-  Drtf(-Olflo»ion 
DrtH-Oiffo»(on(V  .Kfl  lV) 


Physical  Sciences  Research  Laboratory 


Computotfona/ 
-  Materials  Group 


•  Centroid  induced  V,  shift  is  approximateiy  SOmeV  for  the  50nm  MiT  MOSFET 
{t<,ga2nm)  and  lOOmev  for  the  25nm  MiT  MOSFET  (t„,*1.5nm)  as  compared  to  drift 
diffusion  data 

•  increasing  DIBL.  sub-threshold  slope,  and  are  observed  with  decreasing 
(and  tox)  indicating  the  centroid  shift  may  be  compromising  resistance  to  short 
channel  effects 

•  Sub-threshold  slope  and  l^f,  increase  are  also  influenced  by  increasing  sub¬ 
threshold  tunneling  with  decreasing  Lpo,y ...  further  study  into  relative  effects  of 
sub-threshold  tunneling  and  centroid  shift  Is  required 

•  Drift-diffusion  model  show  excellent  agreement  with  calibrated  quantum  transport 
theory  in  predicting  post-threshold  behavior 

Conclusion 

/^proprlately  calibrated  effective  potential  corrections  applied  to  drift  diffusion  models 
will  provide  all  the  essential  physics  required  for  TCAD  device  simulation  for  the 
forseeable  future 


ComputttionMf 

-  Physical  Scisncss  Rsssarch  Laboratory  »  Uatorials  Group 
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Electron-phonon  quantum  kinetics  Electron-phonon  quantum  kinetics  Electron-phonon  quantum  kinetics 

memory  effects  and  Markov  limit  two- particle  correlations  two-particle  correlations 


fourth  Born  approximation,  diagonal  terms: 


Carrier-phonon  Interaction  Carrier-phonon  interaction  Equations  of  motion 

l-b3nd  model  2-band  model  _  distribution  functions  and  polarization 
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Carrier-phonon  interaction  2-band  model  with  constant  field  JntracoHisional  field  effect 

2-band  model  equations  of  motion  2-band  model,  quantum  wire 


spatially  inhomogeneous  system  Spatially  resolved  kinetics  Kinetics  of  charge  separation 

Wigner  representation  density  and  mean  energy  coherent  phonons 


Ch.  Lienau  et  at.,  phys.  stat.  sol.  (a)  178,  471  (2000) 


Carrier  trapping  in  dot  states  Carrier  trapping  in  dot  states  Carrier  trapping  in  dot  states 

_ shallow  smooth  potential _  shallow  smooth  potential  _ resonant  case _ 

initial  energy;  18  meV,  width;  7.5  meV  _ .  .  . — - ,  initial  energy;  18  meV,  width;  7.5  meV 


U13)  ,d  Airausp  1  lAa«u)  jod 


,d  Aiiaiiop  t  (Aaui)  Kni  (^.uo)  ,d  A}iau*p  |  (abui)  ^od  A)l«tiap  )  (a*ui)  -Ktd 


=>  Non-monotonous  behavior  of  occupations 
=?►  Trapping  into  superposition  of  eigenstates 


Carrier  trapping  in  dot  states  Quantum  kinetics  in  semiconductors 

square  well  potential  conclusions 
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Modeling  Ballistic  Current  Flow  in  Carbon  Nanotube  Wires 


Topics  Studied 


M.  R  Anantram 

Nanoelectronics  &  Device  Modeling  Group 
NASA  Ames  Research  Center 

Moffett  Field.  CA 
U.  S.  A. 


Nanotube  >  carbon  nanotube 
HETAl.  METAL 
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*  Role  of  Defects 

*  Bragg  reflection 

*  Metal-nanotube  coupling 

-  wave  vector  conservation 

-  Armchair  versus  Zigzag 


Phonons  /  Mechanical  Deformation 


'  Bragg  reflec.:  Intrinsic  mechanism,  which  exists  even  in  an  ideal  situation 


Curront-carrying  capacity  of  carbon  nanotubes 


•  At  higher  snergiei,  Conductance  =  (<  IkQ) 

n 

An  important  question  is  if  subbands  at  higher  energies  be  accessed 
to  drive  large  currents  through  these  molecular  wires? 


Bias  at  which  electrons  are  injected 

into  non  crossing  subbands  is  aEmc 


7A  56A 


size 

mjyjj 

AE„c(«V) 

1., 

0.98 

0.5 

For  example,  in  a  (20,20)  nanotube  electrons  are  Injected  into  over  20 
subbands  at  an  energy  of2.5eV. 


Experiments  by  various  groups  have  shown  that  the  differentia!  resistance 
of  a  nanotube  decreases  by  small  amounts  with  increase  in  applied 
voltage.  i^.,the  current  carryjnq  capacity  does  not  increase  better  than 
ordinary  resistor  with  applied  bias. 


The  maximum  conductance  if  the  Fermi  energy  Is  at  2.5  eV  is  = 


One  exception  is  the  experiment  by  Frank  et.  al  in  Science  2AC  (1998) 
which  showed  a  rgode^t  decrease  in  differential  resistance  with  bias. 


Main  experimental  features  are: 


'  Vapplied  200mV,  conductance  =  Go 

*  Vapplied  >  200mV,  conductance 
increases  modestly  to  about  1.75  Gq 


For  large  diameter  nanotubes  such  as 
used  in  the  experiments,  the  non- 
crossing  subbands  open  at  about  100 
meV. 

Further  there  are  about  80  subbands  at 
an  energy  of  about  2  eV. 


•  Yet  the  conductance  is  only  ~  1 .75  Gq 


Semiclassical  picture  is  useful 


•  Tunneling  distance  (Xtunnei)  Screening  length 

•  Barrier  height,  2AE|«|q 

•  Scattering  and  Defects 

•  AEnc  a  1/Diameter.  So,  the  importance  of  Zener  tunneling 
increases  with  Increase  in  nanotube  diameter. 


Model 

•  pi-orbital  based  tight  binding  calculation  iPhys.  Rev.  B  SS.  4882  (1998)] 

•  Ideal  contacts  -  reflection  less  contacts 


*  Electrostatic  potential  drop 

•  Linear 

•  Exponential  (Screening  length,  Lgc) 


•  L=2400  A,  L,c=6.  50,  500  A 


L=10  A 

dl/dV  »  4e2/h  for  Vg  <  2AEnc 

Threshold  changes  with 
diameter 

Barrier  height  (AE|4c)  decreases 
with  increase  In  diameter 

Total  Current  increases  with 
increase  in  diameter 

dl/dV>0forVa>  3.1  V. 
except  for  the  (5,5)  nanotube 

(5,5)  nanotube  AE|^c==  1.9eV 


The  differential  conductance  is  NOT  comparable  to  the  increase  in  the 
number  of  subbands. 

For  a  (20,20)  nanotube,  there  are  35  subbands  at  E  =  ±  3.5V. 


•  e-e  and  e-p  scattering  are  not  included 


Coupling  of  carbon  nanotubes  to  metallic  contacts 


Bragg  reflection  severely  limits  the  current  carrying  capacity 


The  crossing  metallic-type  bands  conduct  current. 

The  non  crossing  semiconducotr  type  bands  do  not  conduct  current 


Large  diameter  nanotubes:  non-crossing  bands  will  partially  conduct 
due  to  Zener-type  tunneling. 


The  electronic  properties  of  nanotubes  are  closely  related  to  chirality. 

The  metallic  versus  semicondting  nature  of  nanotubes  and  the  bandgap 
change  with  deformation  depend  on  chirality. 

Questions: 

•  Is  there  a  preferable  nanotube  chirality  to  maximize  current  flow? 

•  Role  of  wave  vector  conservation? 

•  Explain  experimentally  observed  scaling  of  conductance  with  contact 
length 


J 


NANOTUBE 


■f;  •  \  V  electrical  contact  t 
f  7;;'"  length  >7, -77 : 


Parameters  that  influence  current  flow: 

•  Strength  of  coupling  to  metal 

•  Length  of  metal-nanotube  contact 

•  Defects 

•  Metal  Fermi  wave  vector 


How  do  we  model  the  system? 

•  Same  code  as  in  our  MOSFET  work 

•  n  electron  tight  binding  model 

•  Metal  is  modeled  as  a  free  electron  gas  (kp) 


'  TRL  +  TML  +  fLL“2. 

•  Phys.  Rev.  B,  v.58.  p.  4882  (1998)  and  v.  61,  p.  14219  (2000) 

•  Compute  self  energy  due  to:  (i)  metal  &  (li)  semi-infinite  CNT  leads 

Calculation  of  the  total  transmission  from  nanotube  to  metal,  T^l  are 
presented  for  armchair  and  zigzag  nanotubes: 


Scattering  rate 


n  is  an  integer  and  <}>  is  wave  func.  of  atoms  in  a  1D  unit  cell 
The  scattering  rate  from  metal  to  nanotube  (Born  approx.)  is, 

1/T  a  I  <  I  I  «>„,  >  |2 

5(kx-l<x  )  K  I  V|Y,.nt  I  •I’m  ^  1^ 

This  implies  that 

•  kx  is  conserved 

•  ky  conservation  is  relaxed  due  to  finite  width  of  contact  area 

As  a  result  of  the  difference  in  ky  corresponding  to  E=0  (Fermi  energy), 
important  differences  should  arise  as  a  function  of  chirality. _ 


ARMCHAIR 

ZIGZAG 

E=0  at  ky  =  2)i/3ao  =  0,85  A'^ 

Metal  with  kp^rmi  '^0-85  A*^ 
couples  weakly 

E=0  at  ky  =  0 

No  threshold  for  kper^i 

•  Note  that  threshold  value  is  the  value  of  kf  below  which  the  transmission 
does  not  scale  with  contact  length 

The  main  differences  between  the  armchair  and  zigzag  case  are: 

•  Threshold  value  of  kf  is  =0.85  A'^  for  armchair  nanotubes  (see  kf  =  0. 
75A'^ )  and  is  0  for  zigzag  ftanotubes  (see  kf  =  0.4A’^  ). 

•  Beyond  the  threshold  kf,  transmission  increases  with  contact  length  as  seen 
in  experiment  by  Tans  et.  al..  Nature,  vol,  386,  474  (1997) 


•  For  zigzag  tubes.  Tlm  is  small  for  Jty  <  1 .2  A‘^  as  a  result  of  the  large 
angular  momentum  in  the  circumferential  direction. 


Diameter  dependence  of  Transmission 


The  previous  values  of  T  are  small  compared  to  the  maximum  possible 
value  of  2.  Two  possible  scenarios  to  Increase  T  are: 


0.1  0.3  0.5  0.7  0.9 
Diameter  (nm) 


Transmission  decreases  as  the  diameter  increases 
because  the  nanotube  tends  to  a  graphite  sheet. 


*  Large  contact  length  •  Small  coupling 


electrical  contact  jv,  y 
length  '1  ‘  * 


*  Small  contact  length  -  Large  coupling 


electrical  contact 
length 


For  nanoelectronics,  the  second  option  (right)  is  better. 

We  model  a  contact  length  of  30  unit  cells  {72  A  for  armchair  and  125 
A  for  zigzag  nanotubes),  and  vary  the  coupling  strength.  The  main 
results  of  the  calculation  are  rather  surprising  and  are  presented 
below: 
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Coupling  Strength  (eV)  /  Contact  Length 


•  Armchair:Transmission  is  pinned  at  values  I  W 

close  to  unity  for  metal  kf  of  0.9  and  1.2  A'^  \j  ^ 

•  Zigzag:  Transmission  increases  to  two  j  - 

(maximum  possible  value)  A  A 

•  At  small  coupling  strengths,  transmission  m  1^  k 

increases  more  rapidly  in  the  armchair  case  | 

Metallic-zigzag  are  preferable  from  metal-nanotube  coupling  viewpoint 


0.85  A"' 


(|)W+ (fw +  <!>(•)  =  0  ofk. 


:  metal  kf=0,90  A'^  ->  distance  between  nodes  6  A 

;  Gold  kf=1 .20  A*''  distance  between  nodes  3.4  A 

i  Aluminum  kf=1.75  A'^  distance  between  nodes  2  A _ 

Scattering  rate  from  Meta!  to  Nanotube  a  <(tinanotube  I  •^coupling  1  Vmetal> 
<+“1  Hcoupling  I  Vmetai^  ~  0-  Hence,  Transmission  is  pinned  close  to  unity. 

•  Side-contacted:  zigzag  nanotube  are  more  desirable  (curvature) 

•  Larger  metal  Fermi  wave  vector  helps. 


Conclusions 


dl/dV  versus  V  does  not  increase  in  a  manner  commensurate  with  the 
increase  in  number  of  subbands. 


The  increase  in  d!/dV  with  bias  is  much  smaller  than  the  increase  in  the 
number  of  subbands  -  a  consequence  of  bragg  reflection 

Requirement  for  axial  wave  vector  conservation: 

ARMCHAIR  ZIGZAG 

cut-off  ~2iiI33q^  0.85 A  ^  cut-off  Kpgrfjil  -  0 

Our  calculations  show  an  increase  in  transmission  with  lengOi  of  contact,  as  seen 
in  experiments. 

It  is  desirable  for  molecular  electronics  applications  to  have  a  small  contact 
area,  yet  targe  coupling.  In  this  case,  the  circumferential  dependence  of  the 
nanotube  wave  function  dictates: 

•  Transmission  in  armchair  tubes  saturates  around  unity 

•  Transmission  in  zigzag  tubes  saturates  at  two 


WHAT  ARE  CNT? 

STRIP  OF  A  GRAPHENE  SHEET  ROLLED  INTO  A  TUBE 


A  POSSIBLE  LOOPHOLE  IN  THE  THEOREM  OF 

BELL 


Karl  Hess^  and  Walter  Philippi 


^  Beckman  Institute,  Department  of  Electrical  Engineering  and  Department 
of  Physics, University  of  Illinois,  Urbana,  II  61801 
^  Beckman  Institute,  Department  of  Statistics  and  Department  of 
Mathematics,  University  of  Illinois,Urbana,  II  61801 


Abstract 

The  celebrated  inequalities  of  Bell  are  based  on  the  assumption  that  local  hidden 
parameters  exist.  When  combined  with  conflicting  experimental  results  these  inequali¬ 
ties  appear  to  prove  that  local  hidden  parameters  cannot  exist.  This  suggests  to  many 
that  only  instantaneous  action  at  a  distance  can  explain  Einstein,  Podolsky,  Rosen 
(EPR)  type  of  experiments.  The  Bell  inequalities  are  based  on  a  simple  mathematical 
model  of  the  EPR-  experiments.  They  have  no  experimental  confirmation  since  they 
contradict  the  results  of  EPR-experiments.  Furthermore,  in  addition  to  the  assumption 
that  hidden  parameters  exist,  Bell  tacitly  makes  a  variety  of  other  assumptions  which 
could  contribute  to  his  being  able  to  obtain  the  desired  contradiction.  For  instance, 
Bell  assumes  that  the  hidden  parameters  are  governed  by  a  single  probability  measure 
independent  of  the  analyzer  settings.  Moreover,  we  argue  that  the  mathematical  model 
of  Bell  excludes  a  large  set  of  local  hidden  variables  and  a  large  variety  of  probabil¬ 
ity  densities  of  these  local  (in  the  sense  of  Einstein)  variables.  These  exclusions  have 
neither  a  physical  nor  a  mathematical  basis  but  are  based  on  Bell’s  translation  of  the 
concept  of  Einstein  locality  into  the  language  of  probability  theory.  Our  additional 
set  of  local  hidden  variables  includes  time  like  correlated  parameters  and  a  generalized 
probability  density.  We  prove  that  our  extended  space  of  local  hidden  variables  does 
not  permit  Bell-type  proofs  to  go  forward.  This  opens  the  possibility  to  explain  EPR- 
experiments  within  Einstein’s  space-time  continuum,  with  a  finite  velocity  of  light  and 
without  recourse  to  spooky  action  at  a  distance. 


1  Introduction 

Einstein,  Podolsky  and  Rosen  (EPR)  [1]  challenged  Bohr  [2]  and  the  completeness  of  quan¬ 
tum  mechanics  by  designing  a  ’’Gedanken”  experiment  that  suggested  the  existence  of ’’hid¬ 
den  parameters”  and  of  a  theory  that  was  more  complete  than  quantum  mechanics.  The 


1 


EPR  design  was  later  realized  in  various  implementations  [3]  with  experimental  results  close 
to  the  quantum  mechanical  prediction.  These  experimental  results  by  themselves  do  not 
have  any  bearing  on  the  EPR  claim  that  quantum  mechanics  was  incomplete,  nor  on  the 
existence  of  hidden  parameters.  However,  inequalities  derived  by  Bell  [4]  that  are  based  on 
the  assumption  that  local  hidden  parameters  exist,  taken  together  with  the  experimental 
results  that  happen  to  be  inconsistent  with  the  result  of  the  Bell  inequalities,  do  appear  to 
prove  that  local  hidden  parameters  cannot  exist.  This  has  been  discussed  in  great  detail  in 
[5]  and  [6]. 

The  Bell  theorem  is  based  on  a  mathematical  model  of  the  EPR-  experiments.  It  has, 
by  itself,  no  experimental  confirmation  since  its  result  contradict  the  results  of  the  EPR- 
experiments.  The  standing  of  the  Bell  theorem  has  therefore  unique  features  in  the  history 
of  modern  physics:  the  mathematical  model  and  the  theorem  of  Bell  are  taken  to  be  correct 
and  seen  by  many  as  valid  as  the  second  law  of  thermodynamics  is,  while  there  exists  no 
experimental  confirmation.  However,  instead  of  discarding  altogether  a  mathematical  model 
that  contradicts  experiment,  the  contradiction  to  the  experiment  is  used  to  prove  that  the 
basic  assumption  of  the  theorem,  the  existence  of  local  hidden  parameters  is  incorrect.  The 
framework  of  research  that  has  developed  around  the  Bell  theorem  claims  the  necessity  of 
’’gross  non-localities”.  In  simple  words,  the  correlated  spins  of  the  EPR-experiment  are  in 
some  contact  over  arbitrary  space  like  distances  of  our  space-time  continuum  and  if  one  spin  is 
measured  in  one  station,  the  correlated  spin  in  another  station  is  instantaneously  influenced. 
This  contradicts  the  locality  conditions  of  Einstein  and  Einstein’s  very  argument  for  the  lack 
of  completeness  of  quantum  mechanics.  Einstein  called  the  instantaneous  interaction  of  the 
spatially  separated  spins  ’’spukhafte  Fernwirkungen  (spooky  action  at  a  distance)”.  He  did 
not  accept  the  possibility  of  such  spooky  action  and  since  quantum  mechanics  appeared  to 
demand  it,  it  had  to  be  at  least  incomplete.  The  Bell  theorem  and  its  standard  interpretation 
has  turned  the  logic  around.  Its  supporters  now  claim  that  local  hidden  parameters  do 
not  exist  and  cannot  explain  the  EPR-experiments.  Quantum  mechanics  does  agree  with 
experiments  and  spooky  action  at  a  distance  must  be  accepted  as  a  fact  of  nature.  However, 
it  has  been  shown  in  a  serious  of  papers  of  which  we  cite  only  two  of  the  more  recent  [7] 
[8],  that  Bell’s  theorem  does  contain  more  than  self-evident  locality  assumptions.  These 
are  related  to  the  role  of  time  in  the  experiments  and  the  admissibility  of  more  general 
probability  measures. 

We  show  in  this  paper  that  the  assumption  of  the  existence  of  local  hidden  variables  is 
not  the  only  assumption  in  the  proof  of  the  Bell  inequalities.  We  show  that  the  mathemati¬ 
cal  model  of  Bell  excludes  a  large  set  of  hidden  variables  and  a  large  variety  of  probability 
densities  of  these  variables  that  all  fulfill  Einstein’s  locality  conditions  perfectly.  This  ex¬ 
clusion  has  neither  a  physical  nor  a  mathematical  basis  but  is  based  on  Bell’s  mathematical 
interpretation  of  what  Einstein  locality  means  in  terms  of  probability  theory.  Our  additional 
set  of  hidden  variables  or,  as  we  will  call  them,  parameter  random  variables  includes  time 
like  correlated  parameters  and  a  generalized  probability  density  which  is  a  sum  of  what  we 
later  define  as  setting  dependent  subspace  product  measures  (SDSPM’s).  We  demonstrate 
that  Bell  type  proofs  cannot  go  forward  using  our  extended  space  of  hidden  variables. 

The  paper  is  organized  as  follows.  We  first  review  the  theorem  of  Bell.  We  then  analyze 
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the  restrictions  that  Bell’s  proof  puts  on  the  parameter  space  and  probability  measure  and 
show  that  a  much  larger  space  and  a  more  general  probability  measure  can  be  constructed 
without  violation  of  Einstein  locality  conditions.  We  demonstrate  that  a  variety  of  proofs 
of  theorems  similar  to  that  of  Bell  cannot  be  performed  in  this  larger  parameter  space  and 
with  the  more  general  probability  measure,  and  that  these  theorems  and  inequalities  are 
therefore  not  valid  in  this  space.  We  finally  point  toward  a  mathematical  model  that  uses 
this  larger  space  and  permits  the  construction  of  a  hidden  parameter  theory  that  does  agree 
with  EPR-experiments. 

2  The  theorem  of  Bell 

We  first  present  a  short  summary  of  the  work  of  Bell.  In  EPR  experiments,  two  particles 
having  their  spins  in  a  singlet  state  are  emitted  from  a  source  and  are  sent  to  spin  analyzers 
(instruments)  at  two  spatially  separated  stations,  5i  and  52-  The  spin  analyzers  are  described 
by  Bell  using  unit  vectors  a,  b,  etc.  of  three  dimensional  Euclidean  space  and  functions  A  = 
±1  (operating  at  station  Si)  and  R  =  ±1  (operating  at  station  52):  furthermore  A  does 
not  depend  on  the  settings  b  of  station  82^  nor  B  on  the  settings  a  of  station  Si  (Einstein 
separability  or  locality).  Bell  permits  particles  emitted  from  the  source  to  carry  arbitrary 
hidden  parameters  A  of  a  set  Q  that  fully  characterize  the  spins  and  are  ’’attached”  to  the 
particles  with  a  probability  density  p  (we  denote  the  corresponding  probability  measure  by 
p).  Neither  the  parameters  A  nor  the  probability  measure  p  are  permitted  to  depend  on  the 
settings  at  the  stations.  Einstein  separability  is  again  cited  as  the  reason  for  this  restriction. 
The  analyzer  settings  are  changed  rapidly  in  the  experiments  and  do  change  after  emission 
from  the  source.  Therefore  the  source  parameters  and  their  probability  measure  must  not 
depend  on  the  settings  at  the  time  of  measurement.  Bell  further  assumes  that  the  values  of 
the  functions  A  and  B  are  determined  by  the  spin  analyzer  settings  and  by  the  parameters 
such  that: 


A  =  Aa(A)  =  ±1  and  B  =  Bb(A)  -  ±1  (1) 

Thus  Aa(A)  and  Bb(A)  can  be  considered  as  stochastic  processes  on  f],  indexed  by  the 
unit  vectors  a  and  b  respectively.  Quantum  theory  and  experiments  show  that,  for  a  given 
time  of  measurement  for  which  the  settings  are  equal  in  both  stations,  we  have  for  singlet 
state  spins 

Aa(A)  =  -Ba(A)  (2) 

with  probability  one.  Bell  further  defines  the  spin  pair  expectation  value  P(a,  b)  by 

P(a,b)  =  f  A^{X)Bb{X)p{X)dX  =  -[  A^{X)Ay,{X)p{dX)  (3) 

Jo.  J  a 

From  Eqs.(l)-(3),  Bell  derives  his  celebrated  inequality  [4] 

l  +  P(b,c)>|P(a,b)-P(a,c)|  (4) 
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Figure  1;  Light  cone  used  by  Bell.  X  denotes  the  space  coordinate  and  t  the  time 


and  observes  that  this  inequality  is  in  contradiction  with  the  result  of  Quantum  Mechanics: 

P(a,  b)  =  -a  •  b  (5) 

Here  a  •  b  is  the  scalar  product  of  a  and  b. 

The  proof  of  Bell’s  inequality  is  based  on  the  obvious  fact  that  for  x,  y,  z  =  ±1  we  have 

\xz  -  yz\  =  \x  -  y\-l  -  xy  (6) 

Substituting  x  —  >lb(A),  y  ~  =  Aa(A)  and  integrating  with  respect  to  the  measure 

H  one  obtains  Eq.(4)  in  view  of  Eq.(3).  Thus,  from  the  vantage  point  of  mathematics,  the 
Bell  inequality  is  a  straight-forward  consequence  of  the  set  of  hypotheses  and  assumptions 
that  are  imposed. 

3  Extension  of  Bell’s  parameter  space  and  probability 
measure 

We  are  going  to  argue  below  that  Bell’s  parameter  space  is  not  general  enough  and  excludes 
without  any  necessity  a  manifold  of  parameters  that  has  at  least  the  cardinality  of  the 
continuum.  Bell’s  probability  measure  is,  correspondingly,  not  as  general  as  the  physics  of 
relativity  would  permit.  To  show  this  we  start  with  a  discussion  of  the  parameter  space  and 
corresponding  probabilities  out  of  Bell’s  book  [9]. 

Bell  [9]  defines  the  following  parameter  sets  that  are  in  the  backward  light  cone,  as  defined 
by  relativity  and  as  illustrated  in  Fig.  1.  He  lets  N  denote  the  specification  of  all  entities 
that  are  represented  by  parameters  and  belong  to  the  overlap  of  the  backward  light  cones 
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of  both  space-like  separated  stations  Si  and  S2.  In  addition  he  considers  sets  of  parameters 
La  (our  notation)  that  are  in  the  remainder  of  the  backward  light  cone  of  Si  and  Mb  for  S2 
respectively.  Bell  (see  p56  of  ref.  [9])  denotes  the  conditional  probability  that  the  function 

assumes  a  certain  value  with  =  ±1  by  {^a|f^a,  iV}  and  similarly  for  Bb  =  ±1.  Then 
he  claims  that  in  a  local  causal  theory  we  have: 

{Aa|La,iV,Bb}  =  {AalLa.iV}  and  {BblMb,N,A,}  =  {BblMb,N}  (7) 

Eq.(7)  appears  entirely  plausible  as  a  consequence  of  the  finite  speed  of  light:  whatever 
happens  at  station  to  the  result  of  A^  cannot  be  causally  connected  to  the  result  of  Bb 
in  station  S2  within  a  local  theory.  When  the  switching  of  the  settings  is  fast  enough,  the 
probability  that  Ag,  assumes  a  certain  value  must  be  the  same  no  matter  what  value  Bb  might 
assume.  While  this  conclusion  is  undoubtedly  correct  at  a  given  instance  of  measurement, 
Bell’s  use  of  Eq.(7)  as  identical  and  valid  for  all  times  of  measurement  with  a  given  setting 
is  fatally  flawed.  The  reason  is  the  possible  dependence  of  A^  and  Bb  on  time-like  correlated 
parameters  that  may  be  setting  dependent.  The  mathematics  of  Bell- type  proofs  requires 
complete  statistical  independence  of  A^  and  Bb  for  the  whole  set  of  measurements  and 
not  only  at  a  given  time.  It  also  contains  the  assumption  of  identical  La  and  Mb  for  all 
measurements  of  a  run.  This,  however,  cannot  be  guaranteed  because  physical  phenomena 
other  than  the  setting  of  the  polarizer  by  the  experimenter  can  occur  in  the  stations  and 
these  can  be  correlated. 

Consider,  for  example,  two  clocks,  one  in  each  station.  These  clocks  may  have  different 
settings  (e.g.  pendulum  length  and/or  starting  time  etc.).  The  time  that  one  clock  shows 
is  certainly  not  the  causal  reason  for  the  time  of  the  other  clock.  It  is  the  same  physical 
law  that  is  at  work  in  both  stations  and  that  causes  a  correlation  in  the  periodicity  of  the 
processes  in  the  clocks  or  in  some  general  periodic  processes  for  that  matter.  It  is,  of  course 
well  known  that  two  gyroscopes  in  the  two  stations  could  also  be  used  as  clocks  as  they  may 
indicate  the  rotation  of  the  earth.  As  mentioned,  there  may  be  also  other  periodic  processes 
that  cause  correlations  and  these  correlations  may  be  influenced  by  the  settings  a  and  b. 
Although  there  are  clear  analogies  of  gyroscopes  and  spin  properties,  we  do  not  wish  to  push 
this  comparison  too  far.  We  do,  however,  wish  to  point  out  the  dangers  of  using  of  Eq.(7) 
without  proper  caution.  Bell’s  argument  resulting  in  Eq.(7)  does  not  include  the  vital  fact 
that  the  experiments  are  made  in  a  time  sequence  and  that  the  backward  light  cones  change 
and  evolve  with  time.  The  situation  is  illustrated  in  Fig.2  which  shows  that  for  each  instant 
of  measurement  there  is  a  different  light  cone.  Fig.  2  illustrates  our  point  for  measurements 
at  two  different  times  ti  and  t2- 

The  backward  light  cones  contain  sets  of  parameters  La,ti,  Mb,ti,  Nt^  and  La,t2)  Afb^tj, 
ATjj  respectively.  It  is  clear,  that  the  set  La,ti  and  the  set  La,t2  may  contain  setting  dependent 
parameters  A*  with  different  probability  densities.  It  is  also  clear  from  the  discussion  with 
clocks  that  the  sets  La,ti,La,t2,  Mb,t„Mb,t2  etc.  need  not  be  statistically  independent.  This 
has  several  consequences  that  do  not  permit  Bell-type  proofs  to  go  forward.  We  outline 
below  the  most  crucial  problems. 

Bell  uses  combinations  of  Eq.(7)  for  different  settings  in  his  proofs  as  follows  (see  e.g. 
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Figure  2:  Light  cones  at  randomly  chose  times  ti  and  t2  of  the  measurements.  The  clocks 
indicate  time- like  and  setting  dependent  correlations  of  parameters  in  stations  Si  and  S2  (in 
this  example,  the  difference  between  the  times  indicated  by  the  clocks  in  the  two  stations 
stays  constant). 

Bell’s  Eq.(lO)  of  ref.  [9]  on  page  56): 

N,  Bb}  -  {^al^a,  N,  B,]  -  {A^\L^,  N}  (8) 

However,  since  the  measurements  with  setting  b  and  c  are  necessarily  taken  at  different 
times,  this  needs  to  be  written  in  the  form: 

which  is,  in  general,  clearly  incorrect. 

The  above  arguments  also  demonstrate  that  Bell’s  use  of  a  single  probability  density  p(A) 
that  is  valid  for  all  times  of  a  run  of  measurements  is  in  contradiction  to  physical  intuition 
and  facts:  the  parameter  space  related  to  light  cones  changes  and  evolves,  in  general,  with 
time.  To  describe  this  physical  reality  (if  this  word  is  permitted),  one  needs  at  least  to  admit 
a  time  dependence  of  p(A)  i.e.  one  needs  to  replace  p(A)  by 

ti 

In  addition  one  needs  to  include,  again  in  general,  setting  dependent  parameters  denoted 
e.g.  by  A*(t)  in  station  Si  and  by  Ab*(t)  in  station  S2  with 

m  =  Kit)  (11) 

if  b  =  a  in  order  to  make  it  possible  to  fulfill  Eq.(2).  Bell  has  included  into  his  later  proofs 
(after  publication  of  reference  [10])  setting  dependent  parameters.  However,  he  and  everyone 
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else  assumed  that  the  A*  and  A^*  be  statistically  independent.  He  argued  this  independence 
from  the  fact  that  the  parameters  are  in  different  stations  and  he  did  not  consider  time¬ 
like  correlations  as  described  above.  Bell  assumes  that  setting  dependent  parameters  in  the 
stations  analyzers  (called  instruments  by  Bell)  must  be  statistically  independent  or  do  not 
exist.  This  is  explicitly  stated  in  his  book  ([9]  p  38):  ”..it  is  necessary  that  the  equality 
holds  in  (8)  (which  is  equal  to  our  Eq.(l)),  i.e.  for  this  case  the  possibility  of  the  results 
depending  on  hidden  variables  in  the  instruments  can  be  excluded  from  the  beginning 
Of  course,  in  a  run  with  all  different  settings  that  still  would  be  true.  However,  P(a,  b)  is 
evaluated  from  measurements  with  fixed  a  and  b.  Therefore  a  possible  time  dependence  can 
cause  statistical  correlations.  To  visualize  this,  assume  that  the  parameters  A*  and  A^*  are 
identified  with  the  the  hour  pointers  of  two  clocks  in  the  two  stations.  The  clock  in  station 
Si  is  connected  to  a  plane  that  is  perpendicular  to  the  setting  a  and  the  clock  in  station  S2 
to  a  plane  perpendicular  to  b.  Let  the  direction  of  the  pointers  be  idealized  by  unit  vector 
A*  in  station  Si  and  A^*  in  S2  at  each  instant  of  measurement.  Clearly,  these  parameters  will 
exhibit  time  correlations.  Note  that  it  is  of  no  concern  that  the  measurements  are  taken  at 
random  times.  It  is  the  time  correlation  in  the  two  stations  at  any  given  time  that  matters. 

In  addition  to  the  generalization  of  Bell’s  probability  densities  shown  in  Eq.(lO),  one 
needs  a  further  generalization  and  replace  p(A)  by 

(12) 

Of  course,  to  obey  Einstein  locality,  A*  and  A^*  must  be  station  specific  and  can  only  be 
correlated  by  time  like  correlations  i.e.  by  some  relation  to  local  periodic  processes.  It 
is  also  important  to  note  that  the  station  parameters  A*  and  A^*  cannot  be  emulated  by 
the  parameter  pair  a,  A  or  b,  A,  as  always  implied  by  Bell  by  use  of  his  functions  A  and 
B.  The  source  parameters  A  of  Bell  appear  with  a  given  probability  density.  Since  the 
parameters  A*  and  A^*  can  have  different  probability  distributions  for  different  a  and  b, 
that  are  not  related  in  any  way  to  the  parameters  A,  it  becomes  clear  that  the  joint  density 
Pt,(A,  A*,  Ab*)  can  depend,  a  priori,  on  the  setting  vectors  a  and  b.  It  is  irrelevant  that  by 
lucky  coincidence,  the  triples  (a.  A,  A*)  and  (b.  A,  Aj^*)  could  perhaps  be  written  as  a,  A  and 
b,  A  for  some  A  incorporating  A  and  the  station  parameters.  The  probability  density  that 
must  be  considered,  in  general,  for  all  these  parameters  is  therefore  also  different  from  that 
of  Bell  and  must  exhibit  a  time  dependence.  This  implies  the  necessity  of  a  more  general 
probability  measure  that  includes  time-like  correlated  parameters. 

We  have  shown  [8],  that  a  properly  chosen  sum  of  what  we  call  setting  dependent  sub¬ 
space  product  measures  (SDSPM)  does  not  violate  Einstein-separability  and  does  lead  to  the 
quantum  result  of  Eq.(5)  while  still  always  fulfilling  Eq.(2).  By  this  we  mean  the  following. 
The  probability  space  Q  is  partitioned  into  a  finite  number  M  of  subspaces 

Q  =  (13) 

For  given  a  and  b  ,  a  setting  dependent  measure  is  defined  on  each  subspace 

This  measure  can  be  extended  to  the  entire  space  Q  by  setting 

(Mab),„(fli)  =  0  if  j  #  m  (14) 
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The  final  measure  //  is  then  defined  on  the  entire  space  Q  by 

/^  =  1  (15) 

f  ■*m=l 

and  the  index  m  indicates  the  time  correlations.  In  the  above  notation  we  would  have  m  =  ti. 
We  have  shown  in  reference  [8]  that  a  product  measure  can  be  found  such  that 

(Mab)m  ~  (A*a  ^  l^b)m  (1®) 

This  is  of  minor  concern  for  the  main  argument  presented  here  but  does  provide  flexibility  to 
completely  avoid  any  hint  of  spooky  action  within  our  system  of  setting  dependent  subspace 
product  measures  (SDSPM’s). 

It  is  clear  that  Bell’s  proof  does  not  go  through  with  such  a  probability  measure  since 
integrating  Eq.(6)to  obtain  Eq.(4)  works  only  with  a  single  setting  independent  probability 
measure.  In  addition,  one  can  show  in  a  rather  intricate  proof  [8]  that  the  quantum  result 
of  Eq.(5)  can  be  obtained  with  a  probability  measure  as  in  Eq.(15).  In  other  words  hidden 
parameters  are  possible  if  the  parameter  space  is  properly  extended.  We  also  have  shown  in 
[8]  that  the  parameters  that  are  considered  that  way  show  no  trace  of  spooky  action. 

4  Conclusions 

We  have  presented  a  mathematical  framework  that  is  more  extensive  than  that  of  Bell  and 
permits  the  possibility  of  describing  the  spin-pair  correlation  in  EPR-type  experiments  by 
use  of  hidden  parameters.  A  key-element  of  our  approach  is  contained  in  the  introduction 
of  time-like  correlated  parameter  random  variables  that  also  depend  on  the  setting  of  the 
station  in  which  they  influence  the  measurements.  This  leads  in  a  natural  way  to  a  setting 
dependent  probability  measure  composed  of  subspace  product  measures  (SDSPM’s).  Use  of 
such  SDSPM’s  does  not  permit  the  proof  of  Bell  to  go  forward  (nor  any  other  proofs  of  similar 
theorems  known  to  us  as  given  e.g.  in  references  [9]  and  [11]).  We  conclude  that  setting 
and  time  dependent  parameter  random  variables  present  a  possible  loophole  in  theorems  a 
la  Bell. 
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First  principles  modelling  of 
electron  transport  in  atomic-scale  systems 


Antti-Pekka  Jauho 

MIC,  Technical  University  of  Denmark 

1)  Motivation:  Atomic/Molecular-scale  Electronics 

2)  Modelling:  Challenges 

3)  Nonequilibrium  Green  function  formalism 

4)  The  TranSIESTA  ab  initio  transport  program. 

5)  Results  for  atomic  gold  wires. 

6)  A  toy  molecular  device.  . 

Collaboration  of  several  groups: 

Mads  Brandbyge,  Kurt  Stokbro,  Jeremy  Taylor,  MIC,  Denmark 
Pablo  Ordejon,  Jose-Luis  Mozos,  UAB,  Spain 


Molecular-scale  Electronics  -  Why  ?? 

Smallest  functional  unit;  The  single  molecule  (Avirani  &  Ratner,  1974) 

•  “Top-down”  silicon  technology  cannot 
continue  to  shrink  in  component  size 
(dissipation,  tunneling,cost  of  fabs...) 

“End  of  the  road  map  ” 


X 


s 

’  1 


VEAP 


•  “Bottom-Up”:  The  single  molecule  = 
Smallest  functional  unit 

•  Chemists  can  synthesize  organic  molecules 
with  a  range  of  properties 

•  Use  novel  quantum  effects 


I  ■ 


an  asymmetric  moiecvile  {l,4-Bis((2‘-para-mercaptophonyl)- 
ethinyl)-2-acetyl-amino-5-nitro-benzene)  in  between  two  gold 
electrodes. 


Mechanical  Break  Junction  Experiments 


No  Molecules/Solution 


Molecules  in  solution 


“Single”  Molecule  Experiments 


"Conductance  of  a  Molecular  Junction," 
M.A.  Reed  et  al.,  Science,  278, 252 
(1997). 


Conductance  (mS) 


The  Tour  wire:  Mechanism  ??? 


No  explanation  yet 


From  Mark  A.  Reed  and  James  M.  Tou,  Scientific  Am.  June  2000 


Possible  factors: 

A)  Change  in  molecular  states 

B)  Change  in  geometry 

C)  Change  in  charge  states 

D)  Interaction  of  molecular  dipole 
with  external  field 


Only  calculations  for  the 
isolated  molecule 
has  been  published 


Modelling  atomic-scale  conducters:  Challenges 


o-o' 

O’ •0r  o- 

•  Model  a  molecule  coupled  to  bulk  (infinite)  electrodes 

•  Model  the  molecule-electrode  system  from  “first  principles”: 

No  parameters  fitted  to  the  particular  system 

•  Include  finite  bias  voltage/current 

•  Calculate  the  conductance  (quantum  transmission  through  the  molecule) 

•  Determine  geometry:  Relax  the  atomic  positions  to  an  energy  minimum 


Density  Functional  Theory 


Challenge:  Many  interactint 
electrons 


Way  out:  DFT 

I  If  the  electron  density  is  correct  so  is  the  total  energy, 
W.  Kohn,  Nobel  prize  1998. 


m  Ir-r’l 

Interacting  elections 
+  real  potential 


DFT  gives  a  good 

description  of  energy  and  geometry 

of  isolated  molecules  and  perfect 

crystals 


Non-interacting,  fictitious 
particles  +  effective  potentialf^j 
\^f[n(r)],  n(r)  is  the  electron  density 

Schrodingers  equation 
=  H  =  T  +  Ve// 

instead  of  solving  this  we  solve  for  O: 
G(z)  =  (z-H)-1 

We  can  get  the  density  from  G 


Summary  of  important  relations  from 
Nonequilihnim  Green  Fun(;tion  Tlieory 


We  consider  tlie  llainilt  oniau: 

k,(i  k.<x:n 

Current  is  calculated  cis  follows  (Caroli  et  iiL(7()  s),  Meir  ajtd  Wingreen  (1992), 
.)auho  et  al.  (1994)): 


fRojE 


)  4- /Ji(f:)G^(Gfi)]  I 


Here  the  Green  functions  are  matinees  in  the  indices  and  the  linewidth 

functions  T  ^u‘e  defined  as  (here  given  for  t  iine-independent  situation) 

fx^L 


Tinie-i nde[)en( lent  case  (M eir- W i ngrcci i): 

J  -  ^y<kTr{[r^(e}-rV)jG<(q 

-  /)Kur”((0]  [cq.)  -  g"(<:)]} 

-  I  Jdc[jti<)  -  /5i(<)]Ti-{G‘T^'G'r'-}  , 


where  tlie  secoiici  line  liolds,  if  tiie  central  region  can  be  de.s(‘i‘lbe(i  In  terms  of 
an  effective  one-hody  potent uil  (such  fus  in  do'nsiiy-fnnctional  theory). 


Finite  Voltage  Bias/Current 


Challenge:  non-equilibrium 
electron  distribution 


Contact  Regton 


Right  Lead 


Selfconsistent  solution 


G 


r  Greens  Function 


Conductance  =  Transmission 


Conductance  fLandauer-Buttiker^ 

/  =  ^fdEUR{E)  -  fdE))Ttot{E) 
Ttot{E)'.  Total  Quantum  Transmission 
Ttot(E)  =  Tr[ttt](E) 
t(E)  =  (Im  [Eb]  G(B)  (Im  [EJ 


Implementation:  TranSIESTA 


Setup 


Physical  system:  Representation  (PBC  in 


(x.y) 


Example: 

3-atom  Au  chain  connecting 
Au(lOO)  electrodes. 


i  Periodic  Boundary  Conditions  in  (x,y) 
!  3  3  unit  cell,  4  -point  approximation 


CONDUCTANCE 


Benchmark  atomic-scale  conductors: 
Single-atom  gold  contacts 


ps^Good  statistics: 

1  atom=l  conductance  quant 
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DISPLACEMENT  (A) 

•  Brandbyge  et  PRB  52, 8499  (1995) 


Structures  directly  imaged  Single  channel  transport 


OhDishi  et  al.,  Nature  395,  780  (1998) 


Atomic  short  circuit:  Single-atom  contact 


3«atom  Au  chain  connecting 
Au(lOO)  electrodes. 

3<^3  unit  cell,  rf-point  approximation 


-3  -2  -1 

Enengy(eV) 


Expenment: 


DISPLACEMENT  (A) 


Band  structure  of  inflnite  chain 


Transmission  Elgenchannel! 


Benzene- 1,4-dithiolate  coupled  to  gold  chains:  Finite  Bias 


From  He^d  ataL  Science 278. 252 (1997) 


ConducUfWfldiS)  CoftdUCUncediS) 


Summary 

♦  Many  exciting  and  promising  experiments  on 
atomic/molecular  scale  conductors. 

Many  unknown  factors  in  most  experiments. 

V  We  are  developing  a  tool  for  first  principles 
modelling  (TranSlESTA) 


Preliminary  Platinum 


Linear 


8i00 


lO  \  \  li 

1  ^  i  !  i 

—  ■  k  ■  I  ' 

On  !  i 


1  IF' 

I  -  i  tU  •  ^  I  V'; 


-2  -T  0  1 

Energy  (eV)  [E^O] 


Landauer  /  McKelvey 
Approach  to  MOSFET  Modeling 

Mark  Lundstrom 

Electrical  and  Computer  Engineering 
Purdue  University,  West  Lafayette,  IN 

1.  Introduction 

2.  Landauer/McKelvey  Theory  of  the  MOSFET 

3.  The  Ballistic  MOSFET 

4.  Back-scattering  in  MOSFET’s 

5.  Discussion 

6.  Summary 

Lundstrom  additional  information  at:  www.ece.purdue.edu/celab 
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3.  The  Ballistic  MOSFET. 


Ec  (eV) 


Lundstrom 


(Boltzmann  statistics  and  one  subband) 


Purdue 


Lundstrom 


have  a  finite  channel  resistance? 


I  =1^  -  T 

^DS  /D 


3.  The  Ballistic  MOSFET . 


tox  =  1 .5  nm  =  2  x  10’®  cm® 


Vds(V) 


l<a<1.5 

Lundstrom 


100  nm  technolog}/  node..., 

V- ,  -r-5 . -  ;iv<3;"r 


Vdd  =  1V  DIBL=100mV/V 


losion)  /W=qns  Vj. 

Purdue 


3.  The  Ballistic  MOSFET .  comparison  with  measurements . 


Left  =  115/125  nm  technology 


-  ballistic 

(with  measured  Rs) 

Lundstrom 


NMOS;  ~55%  oflimit 
PMOS:  ~  33%  of  limit 

How  can  a  device  with  L^/f 
~  6  mfps  operate  at  ~'50% 
of  the  ballistic  limit? 

Farzin  Assad,  et  al.  (1999  lEDM) 
Dave  Rumsey 
G.  Timp,  J.  Bude,  et  al., 

(1999  lEDM) 

A.  Lochtefeld,  D.  Antoniadis 
(EDL,  Feb  2001) 

Purdue 


Channel  resistanee  vs.  Yes 


•  the  ballistic  l-V  is  readily  computed  and  the  ballistic 
MOSFET  is  readily  understood 

•  present  day  devices  operate  at  ~  50%  of  the 
ballistic  limit 

•  future  devices  will  have  to  operate  even  closer 

•  backscattering  limits  the  performance  of  realistic 
devices 


Lundstrom 


Purdue 


Ec  (eV) 


,  m  .Oil 


4.  Back-Scattering  in  MOSFETs . 


M.S.  Lundstrom,  “Elementary  scattering  theory  of  the  MOSFET,"  EDL,  18 , 361,  1997 

S.  Datta,  et.  al,  “The  MOSFET  from  a  Tranmission  Viewpoint,"  Superlatrtices  and  Microstructures,  1998 

Lundstrom  Purdue 


4.  Back-Scattering  in  MOSFETs Landauer/McKelvey  model 


Landauer  /  McKelvey 
Approach  to  MOSFET  Modeling 

Mark  Lundstrom 

Electrical  and  Computer  Engineering 
Purdue  University,  West  Lafayette,  IN 

1.  Introduction 

2.  Landauer/McKelvey  Theory  of  the  MOSFET 

3.  The  Ballistic  MOSFET 

4.  Back-scattering  in  MOSFET’s 

5.  Discussion 

6.  Summary 

Lundstrom  Purdue 


5.  Discussion.. 


essential  physical  picture . 


velocity 
,  overshoot 


.Ec(x) 


r  =  {l-t)  =  - 


See:  “Essential  physics  of  carrier 
Transport  in  nanoscale  MOSFETs," 
M.  Lundstrom,  etal. 


r^50%=^e^X, 


5.  Discussion.... 


interpreting  simulations 


5.  Discussion,... 

How  to  increase  T? 

. ^  V.',  *-'7  ;  "5  Cf: 

To  reduce  backscattering: 

1)  increase  mfp  l(|i)  " 

^  . 

2)  decrease  i 

\  velocity 
\  overshoot 

\Ec(x) 

[y  . 

Lundstrom 

Purdue 

6.  Summary 

•  Nanoscale  transistor  physics  is  simply  understood  by 
beginning  at  the  ballistic  limit 


•  Scattering  theory  provides  a  simple,  physicai  view  of 
nanotransistors  and  compact,  anaiytical  models 

•  Useful  for  interpreting  simulations  and  guiding  experiments 

•  A  bridge  to  post-CMOS  devices 


Lundstrom 


Purdue 


Diffusion  and  Clustering  of  Impurities  - 
“a  Problem  that  Cannot  be  Ignored” 


Wolfgang  Windl 

Digital  DMA™  Laboratories,  Motorola,  Inc. 
Austin,  TX 


in  collaboration  with 


Benjamin  Liu,  Dejan  Jovanovic,  Mike  Masquelier  (Motorola  Labs) 
Bias  LH^emaga,  Hannes  Jonsson,  Scott  Dunham  (LTW) 


fiverYWhttfe" 


©  2001  Wolfgang  Windl  (Motorola,  Inc.).  All  rights  reserved. 
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MOSFET  Basics 

Metal  Oxide  Semiconductor  Field  Effect  Transistor 


Doping; 

N:  c,  e.g.  As  P:  holes,  e.g.  B 


•  Analog:  Amplification 

•  Digital:  Logic  gates 


(Donor)  (Acceptor) 


Semiconductor  Technology  Scaling 

•  Feature  size  shrinks  on  average  by  12%  p.a.  (f=  0.88) 


•  Chip  size  increases  on  average  by  2.3%  p.a.  {d=  1.023) 


#  of  transistors  /  area  (~  (f  If)  t  by  35%  p.a. 
Switching  speed  (~  1  /y)  T  by  14%  p.a. 


Overall  performance:  t  by  ~55%  p.a.  or 
~  doubling  every  18  months  (“Moore’s  Law”) 


moramaLa 
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MOSFET  Scaling  Challenges 


Solid  Solubility  (cm 


ITRS  Requirements  &  Solubility  Limit 

Solid  Solubility  of  Si  Dopants  ITRS  &  Solubility  Limit 


•Electrically  active  concentrations 
may  be  much  lower 
•“Practical”  solubility  ~  2e20  cm’^ 


How  Are  We  Doing? 


-2 

Sheet  resistance  (H/cm  ) 


•  1  keV  implants  +  optimum  anneals  can  meet  2005  ITRS  requirements  (100  nm). 

•  Further  reduction  in  implant  energy  can  perhaps  meet  70  nm  (2008)  needs, 
but  cannot  meet  needs  beyond  this  node. 


•  Need  to  find  way  to  exceed  2e20  cm'^  concentration  limit  for  n  and  p  dopants,  or 


Ion  Implantation  and  Annealing 


•  Dopants  inserted  by  ion  implantation 

damage 

•  Damage  healed  by  annealing 

•  During  annealing,  dopants  diffuse  fast 
(assisted  by  defects) 

important  to  optimize  anneal 


ht0f!igenc0 


Bverywtief^i* 


Depth  (pm) 


*S.  Srinivasan,  K.  Beardmore,  N.  Jensen 
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TED  &  Deactivation 


What  you  expect:  -- 
Intrinsic  diffusion  | 


Depth  (pm) 


What  you  get: 
TED: 

•fast  diffusion 
•immob.  peak 


intBifIgBfice 


B’vetywhdfB' 


B  Diffusion  in  Bulk  Si:  Qualitative  Picture 


•  B  diffusion  exclusively  mediated  by  Si  self-interstitials* 

•  B  diffusion  limited  by  self-interstitial  diffusion** 


Wandering  I  finds 
an  immobile  B 


^  Short-lived  B/  complex 
moves  a  small  distance 


I  can  be  eaten  by  V\ 


/breaks  away, 
leaving  B  behind. 


♦A,  Ural,  P.  B.  Griffin,  and  J.  D.  Plummer,  J.  Appl.  Phys.  85,  6440  (1999) 

**  W.  Windl,  M.M,  Bunea,  R.  Stumpf,  S.T.  Dunham,  and  M.P.  Masquelier,  Proc.  MSM99 
(Cambridge,  MA,  1999),  p.  369;  MRS  Proc.  568,  91  (1999);  Phys.  Rev.  Lett.  83, 4345  (1999). 


©moTomoLA 
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Reason  for  TED:  Implant  Damage 


Interstitial  assisted  two-step  mechanism: 


Intrinsic  diffusion: 

Create  interstitial,  B  captures  interstitial,  diffuse  together 

Diffusion  barrier:  -^formC-O  ^bind(P-0 

4eV  leV  0.6eV~3.6eV 


After  implant: 

Interstitials  for  “free”  diffusion  barrier  ~  0.6  eV 


W.  Windl,  M.M.  Bunea,  R.  Stumpf,  S.T.  Dunham,  and  M.P.  Masquelier, 

Proc.  MSM99  (Cambridge,  MA,  1999),  p.  369;  MRS  Proc.  568,  91  (1999);  Phys.  Rev.  Lett.  83, 4345  (1999). 
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Reason  for  Deactivation? 

Phenomenological  considerations: 

•  Boron  is  smaller  than  Si  substitutional  B  strains  Si  lattice 

•  Boron  crystal  structure:  Icosahedron  (buckyball),  threefold  coord. 

at  higher  B  concentrations:  New  structures  form, 
bind  &  deactivate  B 

Experimental  findings: 

•  Structures  too  small  to  be  seen  in  EM  only  “few”  atoms 

•  Clustering  dependent  on  B  concentration  and  I  concentration 

formation  ofBJ^  clusters  postulated; 
experimental  estimate:  m  / «  ~  1.5* 

Approach: 

Calculate  clustering  energies  from  first  principles  up  to  “max.”  m,  n 
Build  kinetic  Monte  Carlo  /  continuum  model  from  it 


*S.  Solmi  etal.,  J.  Appl.  Phys.  88,  4547  (2000). 

MOTOmOLA 
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Previous  Work  -  More  Motivation 


De  la  Rubia’s  group  (LLNL) 
predicts  with  ab  intio  based  kMC 
model*  (up  to  B4/2)  a  B 
“activation  window”  for 
annealing  activation 

*M.  J.  Caturla  et  al,  APL  72,  (1998)  p.  2736 

However: 

Experiment*  finds  no  activation 
window,  once  activated  stays 
activated. 

*Mokhberi,  Griffin,  Plummer  (Stanford). 


Boron  active  (%) 


Calculation  of  the  Clustering  Energies 


•  DFT  plane-wave  code  VASP  (Technische  Universitat  Wien), 
LDA  and  GGA 

•  Simulation  cells  of  64  atoms  (converged  ~7%  vs.  216  atoms), 

=  230  eV,  4^  k-points 

•  Relax  many  different,  “guessed”  initial  structures  for  each 
cluster;  dangerous,  can  miss  groundstate! 


Q 
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B-I  Cluster  Structures  from  Ab  Initio  Calculations 


iftt$nrget'ie0 


evwyiv/ier*' 


B12V- 


X.-Y.  Liu,  W.  Wind!,  and  M.  P.  Masquelier, 
APL  77,  2018  (2000);  similar  to  Lenosky  et  al., 
APL77,  1834  (2000). 
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Influence  of  B  Clustering  on  /  Clustering!? 


(e.g.  Kim,  Wilkins,...  PRB  1997) 


New  ring  structure  for  f: 
•Without B:  Ring  ~  0.2  eV  higher 
than  more  compact  f  structures. 
•With  B:  Ring  has  lowest  energies 
B  promotes  ring  growth!? 


mt$iU§mce 


“Silisils” 

(Demkov,  Windl,  Sankey,  PRB  1996) 


moTomauk 
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B-I  Clustering  Reaction  Paths  &  Binding  Energies 


Catch  / 


(GGA  values,  in  eV) 


1.9  * 


.2.4 


Catch  BI 


.2.2  . 


release  I 


mtalttgefice 


Supply  of  I  (vs.  Bi)? 
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DFT  Flavors:  GGA  and  LDA 


For  some  cases,  energetics  differ  considerably: 


Reaction  binding  energy: 


1.5  eV  GGA 
0.9  eV  LDA. 

moro09aLA 
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Calibration  with  SIMS  Measurements 

•  Continuum  model  recalibration  of  ab  initio  numbers  necessary 

•  Genetic  Algorithm: 

-  Start  with  many  random  parameter  sets  within  boundaries 

-  Select  pairs  of  parameter  sets  (“parents”),  biased  by 

-  Mate  parents  (mix  parameters  randomly),  add  child  to  population 

•  Nearly  all  fitted  parameters  between  LDA  and  GGA  values 


B  fraction  f  Activc  B  fraction 


Activation  Results 


•  No  activation  window  in 
agreement  with  expt. 

•  SIMS  calibration  leads  to 
good  prediction  of 
activation  model 
(hopefully)  reasonably 
physical. 

200  400  600  800  1000  1200  •  Assume  only  active 

Temperature  ("^C) 

Expt:  A.  Mokhberi,  P.  B.  Griffin,  J.  D.  Plummer 
fit  from  SIMS 


every 
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Activation  and  Clusters 


•  LDA,  GGA,  fit  do  not  give  consistent  clustering  pattern; 


•  GGA:  Dominated  by  B37“  cluster  (also 

•  LDA:  Dominated  by  B2/3®  cluster  (also  B3r ) 

•  Fit:  Dominated  by  B4/2® 


moTomaLM 
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Uniqueness  of  Clustering  Model!? 


“You  can  fit  an  elephant  with 
21  parameters.” 
(Ulrich  Schroder) 


0  2  4  6  8  10  12  14  16  1B 

D«pth  (nm) 


Energy  filtered  TEM  B  maps  from  (a)  as- 
implanted  (1  keV,  5el5  cm-^)  and  (b)  900  °C 
annealed  sample;  (c)  bright  field  image  of  (b). 
[Wang  et  al,  APL  77,  3586  (2000)]. 


ioteiUganc^ 


Q 


avarywtxata' 
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Application:  Well  Tempered  MOSFET  (MIT) 


intdlHgance 


ayrarywheta' 


B  profile 


Ic-M 
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Cluster  Energy  and  Band  Structure 


DFT  gap 
too  small 
(0.6  eV) 


intedigenc^ 


^verywtiee^" 


Cells  too  small 
(exciton  wave 
function  ~42  A) 


No  defect 
states  «=» 
more  stable 


B12I7 


2- 
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Evidence  of  Stress  Effect  under  Metal  Gate 


IOOOA  TiN/600A  poly-Si  250A  TiN/1 300A  poly-Si  Poly-Si 


00 

CD 

c/5 

CZ) 


tr 

CD 

o 


Experiments  by  Maiti  et  al.:* 

-  Lgff  function  of  TiN  thickness 

-  Stress  estimate  ~ 

Stress  effect 


TiN  thickness 


*B.  Maiti  et  al,  lEDM,  Session  29,  1998. 


int^ltigetfice 


^s/iarywheriB" 


tM.  Laudon,  N.  N.  Carlson,  M.  P.  Masquelier, 
M.  S.  Daw,  and  W.  Windl,  APL78  (2001). 


MOrOmOLA 

digitaldna 


Stress  Dependence  of  Diffusivities 

•  Diffusivity :  D  =  Dq  exp(-  EIkT) 

•  Under  hydrostatic  pressure:  E  -^E  +  pAV 


Perfect  Si  (or  reference  sys.):  Create  defect  (Bg,  I,  V,...)  AV: 


L.  +  AL., 


*M.  S.  Daw,  W.  Windl,  N.  N.  Carlson,  M.  Laudon,  and  M.  P.  Masquelier  (PRB  July  2001). 


intetUgMOi 
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Cell  Volume  and  Charge  in  Si 

S'-atom  Si  cell  64-atom  Si  cell 


•For  all  defects,  cell  sizes,  etc.,  every  additional  electron  adds 
~  15  to  the  volume 


•  Is  this  real  or  an  artifact  of  the  charged  calculation? 
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Electron  Volume  and  Maxwell  Relations 


•  Maxwell  relation: 


dN  dp 


eap)=eao)+Kp 


•  Pressure  dependence  of 
bands  can  be  calculated 
in  good  agreement  with 
experiment. 


TABLE  IV.  First-  and  sccond-ordcr  coefficients  describing  tlie  dependence  of  the  direct  band  gap  at  F 
under  hydrostatic  pressure  =  +  tot  Si,  Gc,  and  GaAs.  The  experimental  results  arc 

from  Goni.  Syasscii,  and  Cardona  fRcf.  1). 


/Ty  a  (meV/GPa)  b  (meVVGPa*) 


Semiconductor 

Thcor. 

Expt. 

Tlieor. 

Expl, 

Thcor. 

Expl. 

Si 

.1.273 

100.8 

0.05 

Gc 

-0.084 

0.795 

125.4 

121 

0.2 

0.2 

CaAs 

0.41 

1.43 

99.1 

108 

-0.1 

-0.1 

Alouani  and  Wills,  PRB  54  (1996). 


•  We  calculate  again  15  for  the 
electron  volume  in  Si 

seems  to  be  real 

•  Ferry*  ~  5  nm  (de  Broglie). 


Pressure  (khar) 
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*D.  K.  Ferry  and  H.  L.Grubin,  IEEE  1998. 
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Conclusions  &  Outlook 

•  Ab-initio  methods  are  being  used  to  provide  predictive  capability  for 
dopant/defect  profile  evolution. 

•  With  recalibration,  boron  clustering  model  can  predict  quantitatively 
activation  and  SIMS  over  a  wide  range  of  temperatures  and  annealing 
times. 


•  Dopant  deactivation  can  easily  be  50%;  concentrations  >  2e20  cm‘^ 
not  achievable.  Deactivation  can  easily  change  by  10s  of  percents. 


•  TED  can  change  junction  depth  substantially  (factor  2  even  at 
“good”  annealing  conditions)  short  channel  effects;  resistance. 


•  We  calculate  the  volume  of  an  electron  in  Si  to  be  15  A^. 
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